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PREFACE 


This thesis entitled “Certain Integral and Fourier Series Equations 
and Their Applications” is an outcome of the researches carried out by me 
since 1997, under the kind supervision of Dr. R.C. Singh Chandel, Reader 
and Head of the Department of Mathematics, D.V. (P.G.) College, Orai. Now 
this work is being submitted for the award of Ph.D. degree in Mathematics. 


The thesis comprises of seven chapters. First chapter is introductory 
giving various definitions and sources of integral and series equations. 
Second chapter deals with recent literature available on the topic. Third 
chapter finds solution of simultaneous dual integral equations. Triple integral 
equations are delt in chapter Four. Chapter Five and Six treat some series 
equations. Finally in chapter Seven we have given application of triple 


integral equations in a crack problem of elasticity. All references are included 
in the end. 
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CHAPTER -1 



INTRODUCTION 

During the last seven decades the special type of integral equations 
have been proved to be very useful tool in the analysis of mixed boundary 
value problems in Mathematical Physics. The same is also true for series 
equations such as dual series equations, triple series equations, etc. 

Integral equations occur in many fields of mechanics and mathematical 
physics. They also represent the formula for the solutions of differential 
equations. A differential equation can be replaced by an integral equation 
with the help of its boundary conditions. The solutions of these integral 
equations satisfy these boundary conditions. 

1.1 SOURCE OF DUAL INTEGRAL EQUATIONS 

One of the classic problems of mathematical physics is that of 
determining the potential in the field due to a circular disc at unit potential 
placed with its plane parallel to and equidistant from two earthed infinite 
parallel plates. The integral form of the solution found is well suited to 
numerical computation when the plates are fairly far apart. If we use 
cylindrical polar co-ordinates (r,z) we many take the disc to be specified by 
z=0, r<l and the plates as z = ±h. The potential has to satisfy Laplace’s 
equation 


1 


( 1 . 1 . 1 ) 



d 2 w 15V 0 2 v 

1 j 

dr 2 r dr dz 2 


= 0 


with, since there is symmetry about the plane z-0. 


V = L 
dV 


dz 


= 0 , 


0<r <1 
r > 1 


when z = 0 


( 1 . 1 . 2 ) 


and 

V = 0, 0<r<oo when z = h (1.1.3) 

The Hankel transform of V denoted by V has to satisfy the ordinary 
differential equation 

= (1.14) 

& 2 P 

The general solution of this equation is 

V = Ae pz + Be -pz ( 115 ) 

and, since from (1.1.3), V = 0 when z = h, we get A = -B e' ph 


Hence 

V = Be _pz - e - P( z ~ 2h ) 0.1.6) 

The inversion formula for Hankel transform is defined as 
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(1.1.7) 


oo 

f(r)= jf(p)p J n (rp) dp 
o 

This inversion formula then gives 

V = JpB[e- pz -e^ z - 2h )]j 0 (rp)dp (1-1-B) 

o 

and substitution in the boundary condition on z = 0, equation (1.1.2) yields 


J pB[l - e 2hp ]j 0 (rp)dp = 1, 0 < r < 1 

o 

Jp 2 b[i + e~ 2hp ] J o (rp)dp = 0, r <1 

o 


(1.1.9) 


This is a set of dual integral equations, where B is the unknown 
function to be determined. 

Thus there are several problems in potential theory which can be 
reduced to dual integral equations. 

1.2 SOURCE OF DUAL SERIES EQUATIONS 

If we make use of series solution of Laplace’s equation we are led to a 
paid of dual relations involving series. For instance, suppose we wish to find 
the axisymmetric solution V(p,z) of Laplace’s equation in the semi-infinite 
cylinder 0 < p < a, z > 0 satisfying the boundary conditions, 

V(p,z)->0 as z-*co 
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V(a,z)= 0, z > 0 

(1.2.1) 

V(p,0)=f(p) 0<p<l 

(1.2.2) 

8V 

u = 0, 1 < p < a 

& z=0 

(1.2.3) 

Then the harmonic function 


00 

V(p,z) = Z X n a n J o( A 'mP)e _A n Z 
=1 

(1.2.4) 

will satisfy the condition (1.2.1) provided that X\, i = 1 to qo are the positive 

zeros of J 0 (X,a). The conditions (1.2.2) and (1.2.3) are then equivalent to the 

paid of relations 


f>;'a n J„(>. n p)=f(p)i ospcl 

11=1 

(1.2.5) 

00 

2> n J o (^p) = 0, 1 <p< a 

n=l 

(1.2.5) 

where [X n ] is the sequence of positive zeros of Jo(Xa). A pair of equations of 

this type is called a pair of dual series equations. Thus there 

are several 

problems in potential theory [192], which can be reduced to 

dual series 


equations. 

1.3 DEFINITIONS OF INTEGRAL EQUATIONS 


An integral equation is an equation in which an unknown function tobe 
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determined appears under one or more integral signs. If the derivatives of 
function are involved in it, then such type of equation is called an integro- 
differential equation. 

An integral equation is called linear if only linear operations are 
performed upon the unknown function, e.g. 

<j)(u) = J K(u, y)g(y)dy ( 1 • 3 • 1 > 

a 

is a linear integral equation and 

<K U ) = |K(u,y)fe(y)} 2 dy {l.'i.l) 

a 

is a non-linear integral equation. 

As equation of the form 

a(u)tt)(u) = F(u)+ A, |K(u,y)j>(y)dy = 0 (1.3.3) 

Q 

is called the linear integral equation where <j>(u) is unknown function, K(u,y) 
is the kernel of the integral equation, F(u) and a(u) all are known functions, X 
is non-zero real or complex number and O is the domain of the auxiliary 
variable y over which the integration extends. 

Linear integral equations involve the integral operator 
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L= |K(u,y)dy 

Q 

having the kernel K(u,y). For any constants Ci and C 2 it satisfies the linearity 
condition 

LMi (y)+c 2 <t> 2 (y)l = (y)]+ c 2 L[<j> 2 (y)] O - 3 - 5 ) 

where 

L [<i>(y)] = J K ( u ,y)4>(y)dy ( 13 - 6 ) 

Q 

1.4 CLASSIFICATION OF INTEGRAL EQUATIONS 

Linear integral equations are classified into following types: 

1.4.1 Fredholm Integral Equation 

An equation having the domain of integration Q fixed is called a 
Fredholm integral equation, e.g. 

b 

a(u)<t>(u)= f(u)+ Xj K(u, y)<j)(y)dy (1-4.1) 

a 

Fredholm integral equations are of two types: 

(i) When a = 0, then equation takes the form 

b 

F(u) = ljK(u,y)<f>(y)dy, a<u<b (1.4.2) 
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This type of equation, which involves the unknown function <j> only 
under the integral sign is called Fredholm integral equation of first kind. 

(ii) When a = 1, then the equation takes the form 

b 

<j>(u) = F(u) + X,Jk(u, y) <|>(y) dy, a<u<b (1-4.3) 

a 

Such type of equations, which involve the unknown function <|> both 
inside as well as outside the integral sign is called Fredholm integral equation 
of second kind. 

(iii) When a = 1 and F(u) = 0, then the equation reduces to 

b 

<j>(u)=^|K(u,y)<j>(y)dy, a<u <b (L4.4) 

a 

This type of equation is known as the homogeneous Fredholm integral 
equation of second kind. 

1.4.2 Volterra Integral Equation 

An equation having variable upper limit of integration is called 
Volterra integral equation, e.g. 

a(u)(j)(u) = F(u)+ Xj K(u, y)<j)(y)dy ( 1 -4 . 5) 

a 

(i) When a = 0, then the equation takes the form 
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a > “oc 


(1.4.6) 


u 

F(u)= A,jK(u,y)<i>(y)dy, 
a 


This type of equation, which involves the unknown function <j) only 
under the integral sign is called volterra’s integral equation of first kind. 

(ii) When a = 1, then the equation takes the form 



<j)(u) = F(u)+ A,J K(u, y)<t»(y)dy (l- 4 - 7 ) 

a 

Such type of equations which involve the unknown function <j) inside as 
well as outside the integral sign is called the Volterra’s integral equation of 
second kind. 

(iii) When a = 1 and f(u) = 0, then the equation reduces to 


<j>(u) = X Jk(u, y ) <j)(y ) dy (1-4-8) 

a 

This type of equation is known as the homogeneous V olterra’ s integral 
equation of second kind. 

1.4.3 Singular Integral Equation 

An integral equation either having the range of integration infinite or 
kernel has singularities within the range of integration is called singular 
integral equation. Such equations occur frequently in mathematical physice 
and process very unusual properties. One of the simplest integral equation is 
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the Able integral equation: 



The solution of this integral equation is given by 


The integral equation (1.4.9) is a special case of the singular integral 


equation 


where h(x) is a strictly monotonically increasing and differentiable function 
in (a,b) and h'(y) * 0 in this interval. It solution is given by 


similarly, the integral equation 


is such in which h(x) is monotonically increasing and a < y < b. This equatio 


has the solution 


sin arc 

d 

K 

dx 


sin olk 

d 

71 

dx 
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1.4.4 NON-LINEAR INTEGRAL EQUATION 


An equation in which the unknown function appears under an integral 
sigh to a power n(n > 1), is known as a non-linear integral equation, e.g. 


1.5 SPECIAL SETS OF INTEGRAL EQUATIONS 


1.5.1 Dual Integral Equations 


The pair of equations 


where G(y) is a function of y alone, is known as weight function \|/(y) is athe 
unknown function to be determined, K(x,y) is the known function known as 
kernel and f(x), g(x) are given functions in the given intervals, are known as 
“dual integral equations”. In practical applications the kernels are usually 


sin an 

d 

% 

dx 


either of Bessel type J v (xy) or of trigonometric type, sing(xy), cos(xy). 

These equations arise frequently in the solution of boundary value 
problems in which the conditions on one boundary is a mixed one. 

1.5.2 Triple Integral Equations 

Triple integral equations are of two kinds. 

(i) Equations of the First kind: 

The set of equations 

00 

JV(y)k(x,y)dy = f(x), 0<x<a (1 .5.3) 

0 

00 

1 G(y)y(y)k(x, y)dy = g(x), a < x < b (1.5.4) 

o 

CO 

jV(y)k(x,y)ly = h(x), b <x <oo (1.5.5) 

o 

where f(x), g(x) and h(x) are prescribed functions in the given interval, are 
known as triple integral equations of the first kind. 

These equations arise in the solution of mixed boundary value 
problems of methematical physics in which different conditions are imposed 
on three parts of one of the boundries. 
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(ii) Equation of the Second kind 


where f(x), g(x) and h(x) are prescribed functions, are known as “triple 


integral equations” of the second kind 


Such types of equations arise in the solution of mixed boundary value 
problems of mathematical physics in which different conditions are imposed 


on three parts of one of the boundaries 


1.5.3 Higher Order Integral Equations 


In the same way, higher order integral equations may be defined. In 
general, we now define “n-tuple integral equations”. 


(i) n-Tuple Integral Equations of the First kind 


The set of “n-tuple integral equations” of first kind are as follows 


(1.5.10) 



00 

J G(y)\|/(y) k (x, y)iy = $ (x), a M < X < 

0 

i = 2, 4, 6, n and a n =oo 

where fi(x), i = 1,2,3, ... . , n are known functions. Quadruple, 5-tuple and 6- 
tuple integral equations are special cases of these equations. 


(ii) n-Tuple Integral Equations of the Second kind 

The set of “n-tuple integral equations” of second kind are as follows: 


j G(yV(y) k (x, y)ty = gi (x), 

o 


i = 1,3,5, 


ai_j < x < 
n-land a 0 = 0 


(1.5.11) 


CO 

J v(yM x , y)iy = gi (x), a M < x < a 4 ( i . 5 . i 2) 

i = 2,4,6, n and a n =oo 

where gi(x), i = 1,2,3, . . . . , n are known functions. Quadruple, 5 -tuple and 6- 
tuple integral equations are special cases of these equations. 

1.6 SPECIAL SETS OF SERIES EQUATIONS 

1.6.1 Dual Series Equations 

The pair of equations 

■ ' 'CO 

2> n G n P(n,y)=f(y\ 0<y<a (1.6.1) 

n=0 


s 


I 

f; 


i 
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, . . V v, « ? || 


EV(-.»)=A) 


a < y < co 


( 1 . 6 . 2 ) 


where G n is a known function of n and some other parameters, P(n,y) is the 
given polynomial of order n and argument y, and sequence (A n ) is to be 
found, are called “dual series equations” involving polynomial P(n,x). 

These equations have been proved to be very useful tool for finding 
solutions to various mixed boundary value problems of elasticity. 


1.6.2 Triple Series Equations 


These equations are of two kinds: 


(i) Triple Series Equations of the First kind 


The following set of equations: 


! 


i>nP(n,y)=f(y), 0 < y < a 


(1.6.3) 


2 A n G » P ( n '>') = sM' a < y <b 


(1.6.4) 


WJ 

£A n P(n,y) = h(yt 


b < y < oo 


(1.6.5) 


are known as “triple series equations of first kind”. In such types of series 
equations f(y), g(y) and h(y) are known functions and the other symbols have 
their usual meaning. 


14 


(ii) Triple Series Equations of the Second kind 


The following set of equations: 


2> n G n P(n,y) = f(y), 


11=0 


0 < y < a 


2>n p (n,y)=g(y), 


n=0 


a <y <b 


£A„G n P(n,y)=h(y)i 


n=0 


b < y <oo 


( 1 . 6 . 6 ) 


(1.6.7) 


( 1 . 6 . 8 ) 


are known as “triple series equations of second kind”. Here, symbols have 
their usual meanings. 

Both kinds of triple series equations have been used for finding the 
solution of problems of potential theory and crack problems of electricity. 


1.6.3 Higher Order Series Equations 

Similarly we have define quadruple, 5-ruple and 6-tuple series 
equations. In general, we now define “n-tuple series equations”. 

(i) n-Tuple Series Equations of First kind 

The set of “n-tuple integral equations” of first kind are as follows: 

00 

£ A n p (n,y) = fi(y), 

n=0 

where i =1,3, 


a i-i < y < a i 
.n-land a 0 =0 


(1.6.9) 
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aw < y < *i 


( 1 . 6 . 10 ) 


Z A n G a p ( n >y) =f i(y)> 

n=0 

where i = 2, 4, n and a n =co 

where the symbols have their usual meanings. 

(ii) n-Tuple Series Equations of Second kind 

The set of ‘‘n-tuple integral equations” of first second are as follows: 


2>nG n P(n,y) = gi (y), a i _ 1 <y<a i 

n=0 

where i = 1,3, n-land a 0 = 0 


( 1 . 6 . 11 ) 


2 A n P(n, y) = gi (y), a iA <y<a { 

n=0 

where! = 2, 4, n and a n =co 


( 1 . 6 . 12 ) 


Quadryple, 5-tuple and 6-tuple series equations are some special cases of 
these equations. 

1.7 APPLICATION OF INTEGRAL AND SERIES EQUATIONS 


Integral and series equations have been extensively used for solving 
mixed boundary value problems of electrostatics, electricity, diffraction 
theory and fluid mechanics. References for such applications can be found in 
text book of Sneddon [192], Sneddon and Lowengrub [194] and 
Muskhelishvili [135], For recent work on Applications some research paper 
cited in the Bibliography of this Thesis may be seen. 

In the present Thesis we have used triple integral equations for finding 
solution of a crick problem in Chapter Seven. 
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CHAPTER -2 


LITERATURE SURVEY 

In this chapter, we give the historical development of the work done in 
solving different types of integral equations and Fourier series equations. 

2.1 INTEGRAL EQUATIONS 

2.1.1 Dual Integral Equations 

(i) Dual Integral Equations Involving Bessel Function 

Earlier research workers actually took interest on dual integral 
equations after the publication of the book ‘An Introduction to the Fourier 
Integrals’ by Titchmarsh [118]. In this text, the following dual integral 
equations involving Bessel functions were considered: 

OO 

Ju a {l + H(u)}\[/(u) J v (xu)du = f (x), 0 < x < 1 (2.1.1) 

0 

oo 

| \|/(u)j v (xu)du = g(x), 1<X <00 (2.1.2) 

o 

with H(u) = g(x) = 0 and a = -1 . 

The range of the solution of the above equations were extended by 
Miss Busbridge [7], The integral equation method for solving the above 
equations was given by Tranter [120], He assumed the suitable integral 
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representation for the unknown function in terms of another unknown 
function and obtained the solution in the form of Fredholm integral equation, 
which can be solved numerically. Gordon [42] obtained the solution of the 
above equations by using certain Sonine’s discontinuous integrals. 

The solution of more general equations was given by Peters [85], He 
took the equations: 

00 

| t a <t>(t) (xt)dt = f(x\ 0<x<l (2.1.3) 

o 

j(t 2 -k 2 ^«j>(t)j v (xt)dt = g(x), 1 < x < oo, k > 0 (2.1.4) 

k 

Noble [72] solved the above equations by applying multiplying factor 
method. Dwivedi [18] obtained the solution of some other equations by 
applying the same method. 

Fan [39] considered some dual integral equations and simultaneous 
equations and solved them by transferring the equations to general system of 
functional equations in the complex domain. He also gave the applications of 
these equations to solid mechanics and fluid mechanics. 

Kalaba and Zagustin [47] solved some dual integral equations by initial 
value problem approach and obtained the solution in the form of a Fredholm 
integral equation. Tanno [117] also considered certain dual integral equations 
as convolution transforms. 
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Srivastav and Parihar [115] applied generalized function concept to 
obtain the solution of certain dual integral equations. Srivastav [116] 
considered dual integral equations with trigonometric kernels and tempered 
distribution. He also introduced L-2 theory to solve trigonometric equations. 

Tranter [119] gave the solution of dual integral equations with Bessel 
function of zeroth order. Sneddon [102] considered dual integral equations 
with Bessel function of zeroth order and reduced them to an Abel-Schlomilch 
type of integral equation whose solution was well known. 

Nguyen Van Ngok and Popov [70] considered the dual integral 
equations connected with Fourier transforms and solved them by applying the 
method of successive approximations. 

Labedev and Uflyand [49] obtained the solution of dual equations of 
zeroth order in the form of Fredholm integral equation with symmetric kernel. 

Recently Chakrabarti [8] considered the three different sets of dual 
integral equations involving Bessel functions of first kind and of order one. 

Aizikovich [4] gave a method to reduce the dual integral equations into 
an infinite algebraic system. The kernels of these equations are the eigen 
functions of a Sturm-Liouville problem for a second order equation in terms 
of a parameter, which is naturally small. 

Cherskii [9] has considered a multidimensional dual equations of 
convolution type 
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(x e n + ) 


(2.1.5) 




u(x)= |k T (x - s)u(s)ds = g(x), 

Q+UQ + 


u(x) = Jk 2 (x- s)u(s)ds = g(x), (x e Q_) 


( 2 . 1 . 6 ) 


where k is defined on the space l(m“ ) of complex valued functions defined 

on M“ with finite norms and M“ (0 < m < n) is the set of points (xi, X 2 , 

, x m , x m+ j, , x„) whose first m co-ordinates are real, while the others 

are integers. Using the factorization method the system is solved by 
quadratures. Moreover, k is considered as an absolutely integrable function, 
vanishing in the exterior of a set W. 

Mandal [58] applied an elementary procedure based on Sonine’s 
integral to reduce dual integral equations with Bessel functions of different 
orders as kernels and an arbitrary weight function, into a Fredholm integral 
equation of the second kind. The result obtained here is more general in 
nature and includes many results, concerning dual integral equations with 
Bessel functions as kernels, known in the literature. 

Considering the three dimensional problems of the theory of elasticity 
with mixed boundary condition with circular lines of the boundary conditions, 
dual integral equations arise. Such types of dual integral equations considered 
by Abramyan [1] are 

00 

JpXp)Jm(i5r>iP = f(4 0<r<a ( 21 - 7 > 

0 
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r >a 


( 2 . 1 . 8 ) 


00 

J#)m(P r )dP = g(r), 

k 

where a = ±1, f(r) and g(r) are known functions in the prescribed intervals. 

\j/(p) is the function to be determined and m = 0, 1, He solved these 

equations by the method of orthogonalisation of equations in the shorter way. 

Rahman [88] found an effective polynomial solution to a class of dual 
integral equations which arise in many mixed boundary value problems in the 
theory of elasticity. The dual integral equations are first transformed into a 
Fredholm integral equation of the second king via an auxiliary function, 
which is next reduced to an infinite system of linear algebraic equations, by 
representing the unknown auxiliary function in the form of an infinite series 
of Jacobi polynomials. The approximate solution of this infinite system of 
equations can be obtained by a suitable truncation. It is shown that the 
unknown function involving the dual integral equations can also be expressed 
in the form of an infinite series of Jacobi polynomials with the same 
expansion coefficient with no numerical integration involved. The main 
advantage of the present approach is that the solution of the dual integral 
equations thus obtained in numerically more stable than that obtained by 
reducing them directly to an infinite system of equations, in so far as the 
expansion coefficients are determined essentially by solving a second kind 
integral equation. | 
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(ii) Dual Integral Equations with Trigonometric Kernels 


The following is the general form of dual integral equations with 
trigonometric kernels: 


fu>(u) S (xu)du = f(x). 


0<x<l 


(2.1.9) 


J \|/(u) Sm (xu)du = g(x) 5 

J me 


( 2 . 1 . 10 ) 


where \j/(u) is unknown function and f(x) and g(x) are prescribed functions. 
Tranter [120] obtained the solution of above equations taking a = ±-^. These 
equations are also solved by Seddon [102] by elementary method. 

Dwivedi [18] obtained the solution of above equations for general 
values of a with various possible combinations of trigonometric functions. 

Singh and Dhaliwal [100] considered the following set of dual integral 
equations with trigonometric kernels. 


|[~ 2^(1 g8)+l e ? -l^Afe)cosfex)d^ = f(x)i 0<x <a (2.1.11) 

J 2E,8 _+ sinh 2^6 V ; V 


J A(^)cos <;x dE, = 0, x > a 


( 2 . 1 . 12 ) 
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(iii) Dual Integral Equations Involving Inverse Mellin Transforms 


Srivastav and Parihar [115] considered the following set of dual 


integral equations involving inverse Mellin transforms 


Erdelyi [36] solved the following set of dual equations 


Cooke obtained the solution of the above equations with f(x) = h(x) = 0 


and later extended the method for f(x) * 0 


Tweed [124] solved the dual integral equations involving the inverse of 
certain Mellin type transforms. Recently Trivedi and Pandey [122] have also 
obtained the solution of dual integral equation involving Naylor’s Mellin type 
transforms by employing Tweed’s method of solution but the equations 
considered by them were quite different from those investigated by Tweed 
himself. Later on Trivedi and Pandey [123] also solved the two pairs of dual 
integral equations involving the inverse of Naylor’s Mellin type transforms. 


(iv) Dual Integral Equations with H-Functions 


Dual integral equations involving H-functions as kernels were first 
introduced by Fox [41], Later Saxena [91] and Saxena and Kumbhat [94] 
extended the result obtained by Fox to more general dual equations. 

Recently Mehra and Ahuja [60] considered the following set of integral 
equations: 

0° r 

r jH(x r ,y r )f(yr) n (dyi)=U(x r ), 0<(x r )<l { 2 . 1 . 11 ) 

o i = l 

where the kernel H(x r ,y r ) is called a multivariable H-function. The second 
integral equation, defined with other boundary conditions, (x r ) > 1, has 
another H-function as kernel with parameters different from that given above 
and involves another function V(x r ) on the right hand side while the prefix (r) 
before the infinite integral sign shows presence of r such operations, the first 
boundary condition 0 < (x r ) < 1 stands for two sided inequalities 0 < xi < 1, 
0 < X 2 <1, , 0 < x r < 1. 

The method used for solution is Fox’s method involving the Laplace’s 
operators L and L' 1 to solve very general integral equation. 

(v) Dual Integral Equations Involving Legendre Function 

Dhaliwal and Singh [17] considered the following pair of dual integral 
equations and obtained a closed form solution. 
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where A(x) is to be determined, P-(i/ 2 ) + ix (cosha) is a Legendre function of 


complex index, f is a real positive constant and f(a) is a prescribed function 


Recently Mandal [59] considered certain dual integral equations 


involving generalized associated Legendre functions of the first kind and 


trigonometric functions as kernels. Solution of these are obtained by using 


properties of the generalized associated Legendre functions and the inversion 


formula for the generalized Mehler-Fok transform involving generalized 


associated Legendre functions of the first kind 


(vi) Dual Integral Equations Involving Hankel Kernel of Order Zero 


Tranter [119] considered the dual integral equations involving Hankel 


kernel of order zero and solved them, which were not derivable by Titchmarsh 


and Busbridge’s solution. Those equations are 



Labedev and Uflayand [49] obtained the solution in the form of 
Fredholm integral equation with symmetric kernel of the equations 


J {l - H(t)}(|)(t) J o (xt)dt = f (x), 0 <x < a (2.1.22) 

o 


j<t>(t)j 0 (xt) = o, x< a (2.1.23) 

o 

Later on Nasim [67] has shown that the dual integral equations with 
Hankel kernel 



GO 

Jt -2a J v (xt)[l + o>(t)j4>(t)dt = f (x), 0 < x < 1 (2.1.24) 

o 

Jt" 2p J ti (xt)i>(t)dt = g(x), x>l (2.1.25) 

o 

where ce> is an arbitrary weight function, can be reduced to a singular 
Fredholm integral equation of the first kind. 

Mehra and Ahuja [60] solved fractional integrals and a Hankel type 
transformation on certain spaces of test functions and of distributions. They 
have shown that the classical relations between fractional integrals and the 
Hankel transformation persist of these spaces, and thus a classical solution of 
certain dual integral equations also persist. 
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(vii) Dual Integral Equations Involving Fourier Transforms 


Nguyen and Popov [70] studied dual integral equations on a system of 
intervals 

i 00 

— f A(£) sin (£x) dS, = f n (x), xel n ,n = l, ,N (2.1.26) 

dxj 

N 

f A(^)cos(^x)l5 = 0, x e R + \ u I n (2.1.27) 

o n = 0 

where A(£) is an unknown function, f„ are given functions and I n = (a n , b„) are 
certain finite non-intersecting intervals contained in R + , where ai < a 2 < ... a n . 

Authors reduced the dual integral equations into an equivalent system j 

of integral equation of the first kind and they gave the criteria for it to have a 

j 

unique solution. This solution can be obtained by the method of successive 
approximations. The authors gave the complete mathematical justification of 
the formal constructions used and constructions in the appropriate function 
spaces. Recently, Nguyen [71] investigated very elaborately the existence and 
uniqueness of the problems for certain dual integral equations involving the 
Fourier transforms of generalized functions. The following dual integral 
equations are considered 

PF" 1 [K(t)u(t)](x) = f (x\ for x = ©, (2. 1 .28) 


P'F -1 [ u(t)](x) = g(x), for x = ©', 


(2.1.29) 
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where o' = R/o, u = F' 1 [u(t)] e s' is a function to be determined, K(t) is the 
non-negative function, f e D'(e o') and g e D'(o') are given distributions on 
o and o' respectively and P, P' are restriction operators to o and o' 
respectively. 

(viii) Dual Integral Equations of Distributional Watson Transforms Type 

Pathak [83] extended the reciprocal formula f = g A h for the Watson 
transform of the function f, defined by g = f n k on the generalized function 
space M' a, b provided the functions k(x) and h(x) satisfy the functional 
equation K(s) H(l-s) = 1, where A denotes the Mellin convolution and K(s) 
and H(s) are Mellin transforms of k(x) and h(x) respectively. Finally he 
applied this theory to obtain solution of a pair of dual integral equations. 

(ix) Dual Integral Equations of Weber-orr Transform Type 

Recently Nasim [68] solved some dual integral equations involving 
Waber-orr transforms 

W~2 k , v b“ 2a p] = f i W a<p<c (2.1.30) 

W v “2 k ,vfe" 2P v(^pl= -h(pl c < p < 00 (2.1.31) 

where k = 1,2,3, ..., v > -1 and vj 7 is an unknown function. Elementary 
methods are used to construct a general solution that contains many known 
results. 


I||}| 
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(x) Dual Integral Equations of Cauchy, Abel and Titchmarsh Type 


Estrada and Kanwal [38] considered the following singular dual 
integral equations and gave distributional solutions. 


Cauchy type 


ctj (x)f(x)+ pj (x)J^^dt + Y! (x)f^dt = g(x\ 


a 2 (x)f (x) + p 2 M J dt - y 2 (x) J dt = g(x\ 


for aj < x < b] 




for a 2 <x <b 2 


(2.1.32) 


(2.1.33) 


where oti(x), a 2 (x), pi(x), p 2 (x), yi(x), y 2 (x) and g(x) are known distributions 
and f(t) is to be determined. 


Abel Type 


ai(x) j^ dt+PlW |-^ dt+ri ( x)& d t 


ft f(t) 

(t-xf 


= g(x), for a t <x <bj 


(2.1.34) 


“ 2 ^^ d, + fc(x) f(^ dt+T2(x) !(^ dt (2.1.35) 


)]A 


v_[(i 
, J ,(t -*y 


= g(x\ for a 2 < x < b' 
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Titchmarsh Type 







up 


J t 2a f(t) J m (tx) dt = A(x), for 0 < x < 1 


(2.1.36) 


| t _2b f(t) J n (tx) dt = B(x), for 1 <x <oo 


(2.1.37) 


where the J n (tx) and J m (tx) are Bessel functions of order n and m. m, n > - !4; 


a and b are real constants; A(x) and B(x) are prescribed functions. 


Recently Eswaran obtained the solution of dual integral equations that 


are classical for diffraction theory 


■ 4 -k 2 A(u)e illx du = f(x), for|x| < 1 


(2.1.38) 


|A(u)e iux du = 0, 


for x > 1 


(2.1.39) 


by reducing them to an infinite system of linear algebraic equations. 


After that more general pair of the form 


Jf v 1 ? v 2 , ,uA(u)e iux du = f(x)t for|x|<l (2.1.40) 


| A(u)e iux du 


forlx > 1 


(2.1.41) 


Jm 

| 


JM 

I 




■ 3 


■¥ 


M 





11 


J : 



where F is some given function with a number of specified properties 


undergoes a similar treatment 


Prabha [87] obtained a formal solution of dual integral equations of 


two variables with two weight functions using the fractional integral operators 


by reducing them to Fredholm integral equations of the second kind 


Veliev and Shestopalov [127] considered the following set of dual 


integral equations and found a general method to solve them 


jh(a)e ±i ^da = 0 : 


These equations are in the class of functions h(ri) for which h(n) 


r|| 1/2 h(rt)eL 2 (k) under certain assumptions on k(rt) and f(r[). They proved the 


unique solvability of this equation 


2.1.2 Simultaneous Dual Integral Equations 


Erdogen and Bahar [37] obtained for the first time the solution of the 
simultaneous dual integral equations involving Bessel functions. 
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These equations are 


00 00 


|Z c ij( x )Vj ( x ) J Hi (xy)dx = Pi (y\ 

On=l 


yell 


(2.1.44) 


00 

JVj( x ) Jix ; (xy)dx = 0, yel 2 (2.1.45) 

o 

where I = 1,2, ...., n. A particular case of above equations for n = 2 is solved 
by Westman [238], Dwivedi [40] and Dwivedi and Singh [49] also considered 
some other sets of simultaneous dual integral equations. 

Narain and Lai [66] considered simultaneous dual integral equations 
involving Maijer’s G-functions of n-variables. 

We have also considered certain simultaneous dual integral equations 
associated with kernel of Fox and simultaneous dual integral equations of 
convolution type in the present thesis. 


2.1.3 Triple Integral Equations 

(i) Triple Integral Equations Involving Bessel Functions 


Tranter [121] was the first person who extended dual integral equations 
to triple integral equations. He considered the equations: 


00 

jV(u)J v (xu)du = 
0 


f f (4 

l h (4 


0<x <a 
b <x <co 


(2.1.46) 
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00 

J u a \|/(u) {l + H(u)}j v (xu) du = g(x), a < x < b (2.1. 47) 

o 

Srivastava [112] considered the following two sets of triple integral 
equations involving Bessel function as kernel. 

I Set 


00 

js A(s)j v (ps)ds = 
0 



0<p<a 

1 < p < CO 


(2.1.48) 


js A(s) J v (ps)ds = -f (p), a < p < 1 


(2.1.49) 


II Set 


j (a,s + ps 2 ) A(s) J v (ps)ds = 


[0, 0 < p < a 

[O, l<p<oo 


(2.1.50) 


00 

|A(s)j v (ps)ds = f(p), a<p<l (2.1.51) 

o 

The triple integral equations are solved by reducing them to ordinary 
differential equation and the solution of this ordinary differential equation 
together with the inversion of Abel integral equation yields Fredholm integral 
equation of the second kind. 

Cooke [14] put the solutions of triple integral equations involving 
Bessel function on sound footing by making heavy use of Erdelyi-Kober 


33 


operators and introduced new operators to obtain the solutions. 

(ii) Triple Integral Equations with Trigonometric Kernels 

Srivastava and Lowengrub [110] solved first time the triple integral 
equations with trigonometric kernel by applying finite Hilbert transform 
technique. Later on Singh [96] considered the different set of equations 


f A(t)cos(xt)dt _ (2.1.52) 

o l 0 ’ x>b 

00 

J A(t) t tanh(xt) cos (xt)dt = p(x), a <x <b (2.1.53) 

o 

and solved them. 

(iii) Triple Integral Equations Involving Inverse Mellin Transforms 

Srivastav and Parihar [115] considered the triple integral equations of 
the type 


| C+lco 

- J.*)p-M. = , ^ p>1 


ffi(px °<P <a 


2m J 

C-100 


(2.1.54) 


C+1CO 


— 7 Js \|/(s)tanap“ s ds = f 2 (p) ? a<p<l 


2m i 

C-100 


(2.1.55) 


They reduced the equations into dual cosine series equations and then 
wrote the solutions by using Tranter’s method. 
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Lowndes [55] extended the operators of fractional integration, 
considered earlier by Cooke, and obtained the solution of the following set of 
triple integral equations. 


M -1 


rfe+s/8) 

r(£+p+s/s) 




x 


0, 0 < x < a 

fi(x), b < x < 00 


(2.1.56) 


M" 1 


r(i+T|-s/g) 
r(l + r[-a-s/a) 



a<x <b 


where a, P, £, r\, 8 > 0, a > 0 are real parameters. 


(2.1.57) j 


Singh [96] considered triple integral equations of inverse Mellin type. 
Recently Tweed [125], considered triple integral equations involving finite 
inverse Mellin transforms. The equations are solved by the method of 
reducing the triple integral equations to a singular integral equation. 

In this thesis we have also considered the two different sets of triple 
integral equations with Legendre function as kernel and solved them using the 
properties of generalized Legendre functions and inversion theorem for 
Mehler - Fock transform. 


(iv) Triple Integral Equations with H-functions 


Saxena and Kumbhat [95] for the first time solved the triple integral 
equations involving H-functions. Recently Mehra and Prabha [61] obtained 
the formal solution of certain triple integral equations containing H-functions 
of n-variables. 
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(v) Triple Integral Equations with Legendre Functions of Imaginary 
Argument As Kernel 

The triple integral equations containing Legendre function as kernel 
were considered by Srivastava [113], Recently Dwivedi considered more 
general triple equations and gave better solution than those given by the 
previous authors. 

(iv) Triple Integral Equations Involving Mehler-Fock Transform 

Srivastava [113] considered the following two sets of triple integral 
equations involving Mehler-Fock transforms. 

I Set: 


J T m(t)A(T)p p (l/2>fiT (coshx)dx = 0, 


0 < x < a 


o L 


J ~ + r A(t) pf^}^ (cos hx) dt = -f(x\ a < x < 1 


(2.1.58) 


(2.1.59) 


I x m(x) A(x)pf (l/2 ^ iT (coshx)dx = 0, 


1 < X < OO 


(2.1.60) 


II Set: 


J t A(t)p P (l/2)+ix (coshx)dx = 0, 


0 < x < a 


(2.1.61) 
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> ' w * ' •’ - ,v ' - 


00 “ | 

J 7 + x2 m ( x ) A(x)p P / 1/2 v fk (coshx)dx = -f(x), a<x<l 


(2.1.62) 


J T A Wp P (l/ 2 )fix (coshx)dx = 0 , 


l<x<«> (2.1.63) 


The above triple integral equations are first reduced to the problem of 
solving an ordinary differential equation. This differential equation is solved 
together with inversion theorem for some variants of Abel integral equation of 
the second kind. At last the problems of solving (triple integral equations in 
reduced to that of solving) Fredholm integral equations of the second kind. 

2.1.4 Simultaneous Triple Integral Equations 

Dwivedi and Sharma [25] have recently considered the simultaneous 
triple integral equations involving H and G- functions of two variables. 
Dwivedi and Singh [24] generalized the problems of H-functions of n- 
variables. Recently Paliwal and Mishra [76] have also considered the formal 
solution of the simultaneous triple integral equations, which are very useful 
for some crack problem in the mathematical theory of elasticity. 

2.1.5 Quadruple Integral Equations 

(i) Quadruple Integral Equations Involving Bessel Functions 

Ahmad [2] has done the extension of triple integral equations to 
quadruple integral equations and those integral equations are as follows: 


; I 

mam 


III 


111 


'111 

✓ 
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Ju' 2a <j>(t) J v (xt)dt = j^ 1 ^ 

Ju' 2P (|)(t)j v (xt)dt = F 2 ^ 

0 IWA 

He solved the above set of integral equations by suing Erdelyi-Kober 
operators and reduced them in the form of Fredholm integral equations of the 
second kind. 

Cooke [15] solved the above equations with general weight function by 
considered some more general operators. Dwivedi, Trivedi and Kushwaha 
[26] generalized the results of Cooke [15], Gupta and Chaturvedi [45], Saxena 
and Sethi [93] have also solved quadruple integral equations with Bessel 

function kernels. 

Recently Prabha [86] considered cretain quadruple integral equations 
involving Bessel functions as kernels and solved them by the application of 
generalized operators of the Hankel transform by Erdelyi-Kober operators of 
two variables. 

(ii) Quadruple Integral Equations Involving Trigonometric Kernel 

Singh and Jain [99] considered quadruple integral equations involving 


0 < x < a 
b < x < c 


a< x <b 
x > c 


(2.1.64) 


(2.1.65) 



trigonometric kernels, which are as follows: 


oo 

Jt A(t)cot(at)dt 
o 


^(a), 0<a<a 

f 3 (a), b<a<c 


( 2 . 1 . 66 ) 
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. .. — • 




J A(t) cos(at) coth(7rt) dt = 


f 2 (a), a < a < b 


c <a<oo 


(2.1.67) 


The above set of integral equations arises in two dimensional steady 
state heat conduction problems. Singh [98] considered integral equations and 
its applications to electrostatics. 


(iii) Quadruple Integral Equations Involving Inverse Mellin Transform 


Recently Dwivedi, Kushwaha and Trivedi [26] have considered the 
following set of quadruple integral equations involving inverse Mellin 
transform: 


m 4 rCi+n-s/or) 1 = jfi (4 0 ^ x < a 

{r(l + ri + a-s/o) ’ J |f 3 (x), b < x <c 


Ffe + s/a) 

r(^+p+s/s) 


82(4 

84(4 


a < x <b 


c < x <00 


( 2 . 1 . 68 ) 


(2.1.69) 


m 


1 


Other authors also considered certain quadruple integral eqautions of 
inverse Mellin transform type (see bibliography). 


(iv) Quadruple Integral Equations with H-Function 


Saxena and Sethi [93] introduced us the Quadruple integral equations 
with H-functions as kernel. Such quadruple integral equations are the 
extensions of the dual and triple integral equations involving same kernel. 


i 
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(v) Quadruple Integral Equations with Legendre Functions of 
Imaginary Argument as Kernal 


Dange and Singh [16] extended the paper of Srivastava [108] for triple 
integral equations with Legendre functions of imaginary argument as kernel 
to the following set of quadruple integral equations 


J A(v) p.(! /2 )+ iv (cos ha) dV = 


fj (a), 0 < a < a 

f s(a), b < a < c 


(2.1.70) 


Jvtanh(7tV)A(v)p.( 1/2 ) +iv (cosha)dV = (2.1.71) 






fm 


J 




c <a <oo 


m 


where A(V) is to be determined. 


2.1.6 5-Tuple Integral Equations 


None of the authors have written a paper on five integral equations 
involving kernel of any kind till date. Dwivedi and Singh [33] have derived 
the solution of n-integral equations from where the solution of five integral 
equations can be obtained as a particular case. 


2.1.7 6-Tuple Integral Equations 


mm 


(i) Six Integral Equations Involving Bessel Functions 


Ahuja [3] considered following set of six integral equations involving 
Bessel functions as kernel 


■ 

1 1 
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44 

0 < x < aj 

44 

a 2 <x < a 3 

.44 

a 4 < x < a 5 

"OQ^ 

to 

"ST 

a 2 < x < a 

< 84(4 

a 3 < x < a 

,86 (4 

a 5 < x <c 






(2.1.72) 


(2.1.73) 


where £ p and £ q are weight functions, ±p * ±q, the values of p and q are taken 
to be 0, 1,-1. The solution of this set is in the form of simultaneous Fredholm 
integral equations. 

(ii) Six Integral Equations Involving Legendre Functions 

Dwivedi, Kushwaha and Gupta [28], obtained the solution of following 
set of six integral equations involving Legendre functions 


co fa (4 0<a<aj 

f<K x ) P?(i/2)+ix( cosha ) dT = \ 44 a 2 < a < a 3 

0 Met), a 4 <cc <a 5 


(2.1.74) 


i « . . Y i 


• 1 p 
ix Jp 

(cosha)a(j>(x)dx 

82(4 

a! <a <a 2 

=< 84(4 

a 3 <a <a 4 

,86 (4 

a 5 <a <a 6 


(2.1.75) 


where f(a), gi(a), i = 1,3,5; j = 2,4,6 are the prescribed functions and <j>(x) is 
an unknown, to be determined. The final solution is in the form of a Fredholm 
integral equation. 


m 


4 
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Dwivedi and Gupta [29] also considered six integral equations with 
generalized Legendre functions as kernel. Those equations are as follows: 


j vWpIJiq^ ( coshx)dx = fj (a\ a H <cc<a j; 


(2.1.76) 


j = 1,3,5 and a 0 = 0 


J\|/(t)rfi-p+iT T “(i-ix x sinh7rx p_( l 1/2 y fix (cosh a)d 


(2.1.77) 


gjCa^au <a <a;,j = l,4,6anda 6 = ■ 


where fj and gj are prescribed functions and V|/(x) is to be determined. The 
final solution is in the form of Fredholm integral equation of the second kind. 


2.1.8 n-Integral Equations 


(i) n-Integral Equations Involving Bessel Functions 


No author has written a paper on n-integral equations involving kernel 
of any kind. Recently Dwivedi and Singh [33] considered n-integral equations 
involving Bessel functions and obtained the solution by reducing them to 
Fredholm integral equations of second kind. 


2.2 SERIES EQUATIONS 


2.2.1 Dual Series Equations 


(i) Dual Series Equations Involving Bessel Functions 




For the first time Cooke and Tranter [12] studied the following set of 
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dual series equations involving Bessel functions 

00 

Z^n 2P A n J v (p^) = F(p), 0<p<l, (2.2.1) 

n=0 


CO 

ZA n J v (pl n ) = 0, 1 < p < a, (2.2.2) 

11=0 

They reduced the above equations to a system of algebraic equations 
which can be easily solved by numerical methods. 

Sneddon and Srivastav [103] considered many other types of series 
equations. 

Srivastav [212] introduced the equations 


A 0 +X>?A n J v fan)=F<pX 0 < p < c, (2.2.3) 

n=l 


00 I 

Z A n Jv(p^n) = G(p) 5 c<p<l (2.2.4) 

n=l 

■ 1 

. ■ 1 

where {A, n } is the sequence of positive roots of the transcendental equation 

■ ' ■ ' ' ' ' ' 

. | 

AJ v (X)+pJ v (A,)= 0 

I 

H and v being the real coefficients with v = -1/2 and p + v > 0. I 

1 

I 

(ii) Dual Series Equations of Jacobi Polynomials 

I 

1 

Noble [72] considered the following set of dual series equations: I 

- : 

; ■ . ! 
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lP«(v,P)A n J„(a,p,x) = F(x)i 


2 A nJ»(«.P»x)=gW 


0 < x < a 


a < x < 1 


(2.2.5) 


( 2 . 2 . 6 ) 




J 

•i fef 
i | 

■ 


sit 


where J„(a, P, x) denotes the Jacobi polynomial and P n (v,3) is a constant 
defined by 


p ( v P ). r (v + n)r(l + q-p + n) 
nl ,P) r(l + a-v + n)r((3 + n) 


(2.2.7) 


He solved the above equations by developing new method of 


multiplying factor. 


Srivastav [114] considered the following set of dual series equations 


y AnPn(q>PX C0S 9) _ p/Q\ n < A < rb 

“^f(a + n + l)f(jP + n + 3/2) ’ - -ty 


( 2 . 2 . 8 ) 


y A n P n (a ? 3Xcose) () 
n= 0 F(P + n + l)f(a + n + 1/ 2) 


4> < e 


(2.2.9) 


where a > -1/2, P — 1 and P n (a,p) (cos 0) denotes the Jacobi polynomials. 
He obtained closed form solution by using Abel integral equation method. 



I HI 


(iii) Trigonometrical Dual Series Equations 


Tranter [120] considered many different kinds of series equations and 
solved them. After some time Tranter [120] gave an exact solution of 
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trigonometric series equations. Noble and whiteman [74] also considered 


other set of trigonometric series equations 


(iv) Dual Series Equations Involving Generalized Laguerre Polynomials 


Srivastava [106] considered following set of dual series equations and 


solved them by applying Sneddon’s [102] method 


(v) Dual Series Equations Involving Legendre Polynomials 


Collins [10] considered the following set of dual series equations 


m+n 


m+n 


where T, 


(cos 0) is the associated Legendre polynomial. He found the 


solution of the above equations in the form of Fredholm Integral equations of 


second kind 


Pathak [82] considered the equations 
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00 A T~ n n n 

£ / “ rPn+p 5 v(x ,-t) = f (x, t), 

n=0 r V H h n + p 

v 2 j 


P n+p,a( x “ t ) = g(x,t), 


0 < x < y (2.2.14) 


n=0 rl n+— +n +p 


y <x << 


(2.2.15) 


where P n (x,t) is a heat polynomial and solution is obtained in closed form. 


Patil [84] studied the dual series equations 


Sr(8 + uKj z » (x ’ k)=fW 0£x<y 


(2.2.16) 


| 0 T^O z » (x ’ k)=gW 


y <x < 


(2.2.17) 


where Z“(x,k) is the Konhouser biorthogonal polynomial. 


(vi) Dual Series Equations with Generalized Bateman K-Functions 


Srivastava [107] considered the following set of dual series equations 


£iw^ K ^ w=fW 0£x<y 


y — .K?W(x)= g(xl y < x <oo 

^,r(2v + a + n + l) 2 ( n+a ) V ’ 5V h y 


(2.2.18) 


(2.2.19) 


where K“(x) is the generalized Bateman K-function. Srivastava [107] also 






§1111 


■if® 
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obtained the solution of particular cases of the above equations. 

2 .2.2 Simultaneous Dual Series Equations 

Dwivedi et al. [31] solved the simultaneous dual series equations 
involving the product of r-Laguerre polynomial. Dwivedi [21] also considered 
the simultaneous dual equations involving H-functions. 

Narain and Lai [65] obtained the solution of simultaneous dual series 
equations involving generalized Bateman K-functions. Recently Dwivedi and 
Shukla [31] considered the simultaneous dual series equations. 

2.2.3 Triple Series Equations §f 


First time triple series equations were solved the Collins [11], He 
solved the following set of equations 

00 

2 (2n + l)C n P n (cos 8) = 0, O<0<a, 3<0 <71 (2.2.20) 

n=0 


00 

2 (l + H n )C n P n (cos0) = f (0) a <0 < 3 (2.2.21) 

n=0 

which are the equations of the first kind and 


5](l + H n )C n P I1 (cos0)=f(0), 0 < 0 < a, 3 < 0 < 7T 

rt=0 


( 2 . 2 . 22 ) 


2](2n + l)C I1 P 11 (cos0)= 0, a<0<3 

n=0 


(2.2.23) 
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are the equations of the second kind. Srivastava [109] considered triple series 
equations involving series of Jacobi polynomials. 

Lowndes [56], Dwivedi and Gupta [30], and Dwivedi, Gupta and 
Shukla [31] obtained the solution of certain triple series equations involving 
Jacobi polynomials. 

Melrose and Tweed [62] obtained the solution of some triple 
trigonometrical series. Parihar [78] obtained the closed form solution of some 
triple trigonometrical series. Lowndes and Srivastava [56] have considered a 
class of triple series equations with Laguerre polynomial kernels and reduced 
them to triple integral equations with Bessel function kernels. The operators 
in the integral equations are for modified Hankel transformations. Hence, 
identities which connect Erdelyi-Kober fractional integral operators with 
Hankel operator are used to reduce the problem to a single integral equation 
involving a modified Hankel operator. Inversion then leads to an exact 
solution, and thus an exact solution to the series equations can be obtained. 

In chapter five of the present thesis we have given a solution of a new 
class of triple series equations with associated Legendre functions which are 
more general than considered as yet some particular cases have also been 
discussed. 

2.2.4 Quadruple Series Equations 

First time Cooke [15] obtained the solution of following quadruple 
series equations: 
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Xa n r» 2 “J v (i n x)= 


GM 


11=0 


a < x <b 
c< x <d 


(2.2.24) 


° <X<S 

n=0 IW), b <X <C 


(2.2.25) 


Dwivedi and Trivedi [22] considered some quadruple series equations 
involving Jacobi polynomials. 

2.2.5 5-Tuple Series Equations 

Dwivedi et al. [32] obtained the solution of various set of five series 
equations involving Jacobi and Laguerre polynomials; and generalized 
Bateman K-functions. 

In the present thesis we have obtained the solution of five series 
equations involving Jacobi polynomials which are still untouched till date. 

2.2.6 7-Tuple Series Equations 

In Chapter Seven of the present thesis we have considered different 
sets of seven series equations involving generalized Bateman K-functions and 
generalized Laguerre polynomials. 

2.3 APPLICATION OF INTEGRAL EQUATIONS IN THE FIELD OF 
ELASTICITY 

It has been seen that integral and series equations are very useful in the 
theory of elasticity, electrostatics, elastostatics, diffraction theory and 
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acoustics. Particularly, these equations are very much useful in finding the 
solution of crack problems fracture mechanics. 

First attempt for theory of cracks was done by Inglis [46] in the first 
quarter of the 20 th century in which he presented the solution of a problem 
within the classical theory of elasticity concerning the equilibrium of an 
infinite body with an isolated cavity. After that Muskhelishvili [64] in 1919 
solved the same problem in much simpler way. 

Griffith presented his first paper [43] in 1920 and second paper [44] in 
1925 in which he rightly considered fundamentals for the theory of cracks. 

Inglis [46] considered the crack in the shape of an elliptical hole. 
Griffith [44] made the minor axis of the ellipse to be zero and reduced the 
crack into a straight line, which is now known as Griffith crack. 

Sneddon [104] developed dual integral equations and applied first 
transform techniques in solving the variety of mixed boundary value problems 
of the potential theory. Muskelishvili [64] developed a method based on the 
theory of Cauchy’s integral for solving the mixed boundary value problems. 
These methods have played an important role in the development of potential 
theory and mathematical theory of elasticity. Sneddon [104] has dealt with 
different methods to solve mixed boundary value problems. 







Recently Tweed and Melrose [126] also considered ‘The out of plane 
shear problem for an infinite sheet with a staggered array of pairs of cracks. 
They used a triple series technique to find closed form expressions for the 


SO 


mode III stress intensity factors of a staggered array of pairs of cracks in an 
infinite elastic solid. 

A 2-dimensional problem of an anisotropic elastic strip having an 
infinite row of Griffith cracks was considered by Misra and Misra [63] by 
using dual integral equations approach. The problem analyzing analytically 
the stress intensity factor, the critical pressure and the energy required to 
open the crack, were studied. Numerical results were also derived. 

Griffith’s 2-dimensional theory of rupture was extended to 
3-dimensional by Sack [90], He considered a disk shaped crack in the interior 
of infinite elastic solid. It is usually known as “penny shaped”. 

Lai and Jain [51] studied the problem of finding the distribution of 
stress and the displacement jn an elastic half plane containing an external line 
crack perpendicular to the free surface of the plane. This problem has been 
formulated in the form of dual integral equations. Finally, the expressions for 
the stress intensity factor and the crack energy are derived numerically. 

Palaiya and Majumdar [75] solved the problem of 2-diffraction of 
harmonic antiplane shear waves by a pair of coplaner parallel rigid strips 
located at the interfaces by using Fourier cosine transforms. The problem was 
reduced to triple integral equations, whose solution was obtained by usual 
methods. 

Srivastava and Dhawan [111] considered the problem of stress 
distrubition due to a Griffith crack at the interface of an elastic layer bonded 
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to half plane. They reduced this problem to a system of simultaneous dual 
integral equations involving trigonometric kernels. Some further crack 
problems have been solved by using simultaneous dual integral equations. 

Lowengrub and Sneddon [54] considered the problem of stress 
distribution in a dissimilar medium when a penny-shaped crack is situated 
along the bounding plane. They converted this problem to simultaneous dual 
integral equations which were later solved by known methods. 

Certain triple integral equations have been found very useful in the 
solution of three part boundary value problems. Lowengrub and Srivastava 
[52] considered certain problems where such equations arise. Several other 
research workers used triple integral equations to solve different crack 
problems. 

Lowengrub [53] considered the problem of stress distribution for the 
two bounded dissimilar elastic half planes containing a pair of coplaner 
cracks at the interface line by using Fourier transforms. He reduced the 
problems to simultaneous triple integral equations and later to a Hilbert - 
problem. 

Roy and Chatter] ee [89] considered the effect of the free surface of the 
stress distribution of an elliptic crack aligned parallel to the free boundary 
and at depth h-below it. The title problem is posed to dual integral equations 
in Cartesian co-ordinate system. By suitable transformation the dual integral 
equations are first reduced to an infinite system of dual integral in cylindrical 
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coordinates. Then they are further reduced by a recently developed technique 
to an infinite system of Fredholm integral equations of the second kind. 

The problem of determining the stress intensity factors in a semi- 
infinite orthotropic elastic medium containing two coplaner cracks parallel to 
the boundary was considered by Mahapatra and Parti [57], The above mixed 
boundary value problem was reduced to triple integral equations which were 
solved by using infinite Hilbert-transform technique. 

Quadruple integral equations have also been used to find the solution 
of certain crack problems. Singh [97] considered the problem of determining 
the stress distribution in the vicinity of a Griffith crack in an infinite elastic 
solid, when the crack is opened by two symmetrical rigid inclusions. 


Recently, Lai [50] considered the problem of elastic half-space under 
torsion by a flat annular rigid stamp in the linear micropolar elasticity. The 
problem is reduced to system of four Fredholm integral equations. 


2.4 APPLICATIONS OF SERIES EQUATIONS IN CRACK 
PROBLEMS OF ELASTICITY 


A large number of problems are encountered in the theory of crack 
problems of elasticity. Some of these problems are being mentioned below 
where series equations have been found useful. 


Parihar and Garg [79] derived the stress and displacement distributions 
when infinite Griffith cracks are presented along the bound line of two 
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dissimilar elastic half planes. The problem is first reduced to the dual series 
equations and then to a Hilbert problem. 

Parihar and Kushwaha [80] considered the problem of distribution of 
stress in an infinite strip containing two Griffith cracks under the action of 
body forces. They reduced the problem to a set of triple series equations and 
obtained the closed form solution. 

Parihar and Kushwaha [81] also considered a crack problem of a strip 
containing Baranbalatt crack. They reduced the problem to dual series 
equations. By symmetry conditions the problem is equivalent to that of an 
infinite row of collinear Baranbalatt cracks in an infinite elastic solid. 

The problem of stress distribution in the presence of infinite row of 

interface cracks, located symmetrically on the central line in a composite strip 

\ 

was considered by Parihar and Garg [19]. This problem is first converted to 
simultaneous dual series equations and finally to Hilbert problem. 

A number of crack problems in the classical theory of elasticity have 
been solved with the help of integral and series equations. A complete 
bibliography of the work done till 1967 can be found in a monograph by 
Lowengrub and Sneddon [54], 

In the present thesis, we have used triple integral equations to solve a 
crack problem of elasticity. 
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CHAPTER -3 

SIMULTANEOUS DUAL INTEGRAL EQUATIONS 

In this chapter we shall discuss solutions of certain simultaneous dual 
integral equations. In sections 3.1 and 3.2, the simultaneous dual equations 
involving H-functions and convolution type have been considered 
respectively. 

3.1 CERTAIN SIMULTANEOUS DUAL INTEGRAL EQUATIONS 

ASSOCIATED WITH KERNEL OF FOX 

In this section the formal solution of certain simultaneous dual integral 
equations involving H-functions is obtained by the method of fractional 
integration. By the application of fractional integration operators, the given 
simultaneous equations are transformed into two others having common 
kernel and the problem then reduces to that of solving one integral equation. 
Since the common kernel obtained is a symmetrical Fourier kernel given 
earlier by Fox the solution then follows easily. 

3.1.1 Introduction 

Fox [41] and Saxena [91], [92], [94] have obtained formal solution of 
certain dual integral equations involving H-functions. In this section we have 
discussed certain simultaneous dual integral equations of more general nature 
associated with the kernel of Fox [41], By the application of fractional 
integration operators the given integral equations are transformed into two 
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others with common kernel and problem then reduced to that of solving one 
integral equation. Since the common kernel obtained comes out to be a 
symmetrical Fourier kernel investigated earlier by Fox, the formal solution is 
readily obtained. 


We define the H-function by means of Mellin Barnes integral [91], in 
the form 


H 


m,n 

p»q 


|(a p ,A p 

'M 



i(ai,A 1 ).... 

•”( a P’A p T 

(bi,Bj).... 

•■(bp.B,) 



n 

71 

j=l 


P 

71 

j=m+l 


r(b ; + b :s) tc r(l - a r - a r s) 

X“ s ds 

r(l-bj -B:s) 7i r(a r +a r s) 

r=n+l 


(3.1) 


where x is not equal to zero and an empty product is to be interpreted as 
unity; p,q,m and n are integers satisfying 1 < m < q, 0 < n <1; A and Bj are 

positve numbers and a and bj are complex numbers (where r = 1, 2, p and j 

= 1, 2, , q) such that no pole of F(bj + Bjs) (j = 1,2,3 m) coincides with 

any pole of T(1 - B r - sA r ) (r = 1, 2, 3, n) i.e. A r (b r + v) * Bj(a r - T] - 1) for 

v, r[ = 0,1,2,.... 


The above definition possesses one advantage that for a suitable 
contour L, and coefficient of x' 5 in the integrand of (3.1) is evidently the 

Mellin transform of H“q n (x). 


An account of the covergence conditions of the integrand (2.1) could 
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be found in the work ofBraaksma [5] and Erdelyi, et al. [34], 


Fox has shown that the function 



TnrH,m 

n 2m,2n 


I O' a m>^mX a m A m ,A m ) 
(b n 3 n ),0-b n -B n ,B n )_ 


= ^ m m 2n (x) 


2m^ m ’ n ^ S ) X ^ S ’ 


(3.2) 


where 


tmAS> wrfe+Bi-sBijSlfaj-Aj+sAj) 


(3.3) 


behaves as a symmetrical Fourier kernel from (3.2) it follows that Mellin 
transform of H^ 2n (x) is Xm,n(s). 

3-1.2 Simultaneous Dual Integral Equations 

The solution of the following simultaneous dual integral equations will 
be obtained here: 


Tjn,m 

11 2m,2n 


XU 


(l a m?^m)( a rn 

(bn.BnMl-Pn-Bn.Bj, 


S a hk f h( u )du 

Ih=l 


= <t>k(x). 


0 <x < 1 


(3.4) 


and 
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IPn^nMl t>n ®n?®n) h=l 


Vk(x), 


X > 1 


(3.5) 


where k = 1, 2, , n; <|>k(x) and \j/k(x) are given, fh(u) are to be determined 

and ahk, bhk are known constants. The H-functions used here are as defined in 
(3.2). In above equations (3.4) and (3.5) the unknown function fh(u) must be 
bounded and integrable. 


If we put M[f(u)] = F(s) and apply Fox’s result [91] to (3.4) and (3.5), 
we obtain 




!■ 


1 r n 

l • J Xm,n (?)x 0 ~ s )ds = <j>k (x), 


0 <x < 1 


(3.6) 


■fl* 

life 

on 

Mi IhJ: 




l n 

rr J X m ,n (s)x -s 2 bhkFh (l - s)ds - y k (x), * > 1 


(3.7) 




where k = 1,2, n and 


XmA S ) =% . 


,r(bi +sBj) ” r(af-sAj) 


rtef+Bi-sBilj-irlaj-Aj+sAj, 


(3.8) 




lm 
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3.1.3 Reduction of (3.6) and (3.7) to Equations with A Common Kernel 


The well known beta function can be expressed as 


> ct 0 + it on the line a = o, 


These conditions may no longer be necessary and the second may be 
relaxed. Replacing x by v in (3.6), multiplying by v c m a m -l( x c m . v c m ) (a k m -a m -l) 


and integrating (which is justified under the assumptions already made) 
through the integral sign with respect to v from 0 to x, 0 < x < 1, we find that 


where %^ n (s) is obtainable 


r ^ 




The operator of fractional integration denoted by R is used in the form 


R[M : m : w (x)] = _^_ x _8+mX+m_1 

r(w) 


-f(x m -v m ) v 8 w(v)d\ 


(3.12) 


The case m = 1 is due to Kober [48] and the result for more general 
case m > 0 has been given by Erdelyi [34], Fox [41] has shown that there is 
no essential difference between the two cases, since 


R[y,§ : m : w(x)] = Ry,m 1 (g + 1)- 1 : lW(x)J. 


(3.13) 


where x m = X, v m = V and w(x) = W(x). 


If in addition, the operator R(x) exists provided w(x)e L p (0,qo), p > 1, y 
> 0, 8 > 1/q. If w(x) can be differentiated sufficiently often, then the operator 
R exists for negative as well as positive values of y. 


For brevity, we write 


R [( a i k _a i)( a i A 1 rl “ij-Aj 1 : w(x)J= Ri[w(x)]. 

R|bj -3f)(p k -Bj)B - 1 -1,B - 1 : w(x)J= Rj[w(x)|. 


(3.14) 


(3.15) 


From (3.14) it is seen that right hand side of (3.11) is equal to 
Rm[<k(x)] with 0 < x < 1. On transforming the equation (3.6) step by step by 




m . 


means of the operators Ri and Rj for I = m, m-1 2, 1 and j = n, n-1. 


u ! 
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We now proceed to transform the integral equations (3.7). Again from 
beta function formula, we have 


where d n - B n ,P n > b n and d n a 0 +b n > 0. In (3.7) let us replace x by v and 


multiply by v 


and integrate with respect to v from 


x to oo through the integral sign. Equation (3.7) takes the form 


(s) is obtainable from 


where x > 1, d 


replacement of f3 


by b„. The second operator of fractional integration 


denoted by K required here is 
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K[y,S : n : w(x)]= ~,x s j(v" - x n ) T 'v“ 6 ~”' +n “ 1 w(v)dv 


(3.19) 


If w(x) L p (0,qo), p > 1 and w(x) can be differentiated sufficiently often 
K exists, provided n > 0, 8 > 1/p, where y can have any positive or negative 
values. 


-bi)3iB i 1 ;B i w(x)j= kj[w(x)] 


(3.20) 


K [( a j -a f)( a ^A j )A.J 1 :w(x)j=k*[w(x)] 


(3.21) 


from (3.20) it is evident that the right hand side of (3.18) is K n [\j/ k (x)], where 


x> 1. 


The successive application of the operators, Kj for I = n, n-1, 2, 1 

and Kj of j = m, m-1, 2, 1 to (3.7) yields the integral equation 


n 

T Jx m,n ( s ) x _S Z a hk F h 0 - s)ds 


2c hk K;[K;....<K 1 K 2 XjvkWl- 


(3.22) 


where Chk are the elements of the matrix [ahk] [bhk]" 


If we set 
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R 1 ^2 •• R n R l R 2 R m[<f>k( x )] 1 


0 < x < 1 


(3.23) 


- 2 C hk K l [k 2 . ...X n R l K 2 R m[Vk( x )l ] X>1 


h=l 


(3.22) 


k = 1,2, ,n. 

Equations (3.16) and (3.22) can be put into the compact form as 


lr n 

-rJXm ; n(s)x ^a^Fj^l-s^s 


2rci 


h=l 


M M [ H 2m,2n( u )] x S X a hk F h(l " s)ds = g k (x) 


27ci 


(3.24) 


h=l 


where M denotes the Mellin transforms. Since is a symmetrical 

Fourier kernel, the formal solution can be readily obtained as 


n 00 


f (v)— 'V'rj,, f T-f n ’ m (inA 1 

th 2. d hk J H 2m;2n l x uj( bn B a ) (l _ bn _ Bn;Bn) 


k=l 0 


g k (u)du = 



k=l o 


(3.25) 


r:|r‘....r.>,r, R,k(xi) J 


(3.23) 


JhS 


n 5 m 

2m,2n 




Bn )>(l~b n -B n ,B n ) 
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Z C hk K l ^ 2 - • K m K iK 2 K m [\|/ k (x)J ] 


h=l 


(3.22) 


where h = 1,2, ,n. 


3.1.4 Interesting Particular Cases 


(1) If we set m = 0, n = 1, B, = B 2 = 1, b x = X k + y k /2, Jjf =y\-\ and use 

the identity [] 


G‘f 2 (x/a,b)=x a+b ' 2 j a _ b (2x) 1 ' 2 


(3.26) 


We find the formal solution of the simultaneous dual integral equations 
j(ux)A. k v k [2(xu) 1/2 ]^a hk f h (u) = (j) k (x), 0<x<l (3.27) 


h=l 


J (ux) kk v k [2(xu) 1/2 JjT a^fjj (u) = \|/ k 


X < 1 


(3.28) 


h=l 


k= 1,2, ,n; 


is given by 


f h (x) = X dj* J. Jv k + 2X, k [ 2 M 1 n ^.[2X k ,v k / 2 - A, k ,l] (3 ,28a) 


h=i 0 


<t> k (u)du + J^]J k k + 2X 


2(xu) 1/2 tchkK[-2^ k A k + v k /2,l,] 


h=l 


V|/ k (u)du, h = 1,2, ,n. 


(3.29) 
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(ii) On the other hand if we put m = 1, n = 0, Ai = A 2 = 1, 


af-1 p k +H k /2,a 1 -l + p k +ji k /2 and apply the transformation [2] 


. k , ,,k 


g m Wb:l=Gsrk i ]; _b ’ 


(3.30) 


g o° 2 ff x /a M= x^ a+b ) /2 J a _ b (2x) 1/2 , 


(3.31) 


It is found that the formal solution of the simultaneous dual integral 


equations. 


J ( ux ) k 1 J v k[2(xu) 1/2 a hkfh( u )du = <j) k (x), 


0 <x < 1 (3.32) 


J( ux ) k 1 J v k[2(xu) 1/2 ]^bhkfh(u)du = \|/ k (x). 


(3.32) 


is given by 


n 1 

f h( x )=2X k J j v k +2k( xu ) _ 

h=l 0 


n _ 0 k k , 1£_ 

** Q y q I 


,l,<j> n (u)du + j(xu) ! J k +26 1 


[2(xu)- 1/2 ]£ c hk K[2q k ,2p k - q k ,1, ¥k (ujduj 


(3.34) 


where h = 1,2, 
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3.2 SIMULTANEOUS DUAL INTEGRAL EQUATIONS OF 
CONVOLUTION TYPE 

This section deals with certain simultaneous dual integral equations of 
convolution type. The solution has been obtained by following the method of 
Tanno [117]. 


3.2.1 Introduction 


Fox [41] and Saxena [91] have obtained the formal solution of certain 
dual integral equations involving H-function as kernel. Tanno [117] obtained 
the solution of more general class of dual integral equations. This class 
involves the case of Fox and Saxena as special cases. 


Erdogan and Bahar [37] and Westmann [128] obtained the formal 
solution of simultaneous dual integral equations involving Bessel functions as 
kernel. Recently Dwivedi [21] has obtained the solution of certain 
simultaneous dual integral equations: 




tEW 1 ' 1 ’: - ' " 

iv : 

Stf w ... 


J n P+n#n 


. (u klu — (j>u (x j 


f 

J n p+I 


Tb hk f h (u)du = \|/ k (x), 


1 







where Hjf (x) is used to denote the function H^ n (Ha i A i > hk zb j Bj and b h k 
are known constant <J»k(x) and \|/k(x) are given functions and fh(x) to be 


determined 


The H-functions of order n used above are of the form 


where ai > 0, ai, Pi are all real, i = 1, 2 


n, the contour C along which 


the integral of (3.39) is taken as the straight line parallel to the imaginary axis 


in the complex s-plane and lie to the left to the line a = <j 0 > -a 


The integral (3.39), taken along the line, coverages if 


and coverages absolutely if 
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i=l i=l 

Our object there is to obtain the solution of more general class of 
simultaneous dual integral equations of convolution type by following the 
method of Tanno [117], 

From the point of view the simultaneous dual integral equations (3.35) 
(3.36) and (3.37), (3.38) can be reduced to the dual convolution transforms 
after exponential change of variables we consider the followings: 

CO n 

j G k(x-t)2 a i lk <i> h (t)dt = f k (x), x > X (3.40) 

-oo h=l 



|H k (x-t)^b hk \|/ h (t)dt = gi c (x), x > X (3.41) 

—oo h=4 

k = 1, 2, 3, , n. 

where X, ahk and bhk are constants and Gk(t) and H k (t) are generated by 
certain meromorphic functions. The explicit forms of these functions are 
defined in the following section. 

3.2.2 Definition of Kernel 

Let us define admissible entire functions Ej(s) (i = 1, 2, 3, 4): 


E i (s)=exp(8 i s)7t 

j=i 




a bj ) 


(s) 

exp—, 

a ij 


(3.42) 


68 


! 




lEki 


where 81, ay are real numbers such that jri — l <oo (i = 1, 2, 3, 4, 

k=l [ a ij J 

j = l ,2, ) and ay < 0, a 2>j > 0, a 3)j < 0, a 4 ,j > 0 (j = 1, 2, ). 


Let (a, b) be an interval included in the corresponding intervals I; (I = 
1 5 2, 3, 4) which contain no zero of Ej(s) except perhaps at an end point. 


Further, if 




00 hck 


E|(a + ix) 
E^a + ix) 


dx <oo for every o in a < a < b , 


7 E 4 (a + ix) 4 

—rz dx < oo for every a m a < a < b, 

L E^ + h) 


lim Eo(a + ix) 

— r \ dT = 0 uniformaly in a < a < b , 

x -» oo E } (a + ix) 










lim 5l|2±i!) = o uniformaly in a <a<b, 
x— >-ooE 3 (cy+ix) 


K = 1,2,3,. 


Jj|;5 1 ; | 

£ §i|ff 


then there exist two function sFk(t) and Hk(t) such that 


, ■ 1 ts® , 

: f '• 


E|(s). 


E,(s) 


Je st G k (t)dt, a<o<b. 


(3.43) 


E 4(s) 

E|(s) 


= Je ^H^t^t, a <o <b. 


(3.44) 


K = 1, 2, 3, , n. 


«! 



where b\ and E| are obtained from (3.42) replacing a 2 , ja 3 , j by a|, ja|, j 

respectively and the intervals converging absolutely. The classical inversion 
formula [] gives F k (t) and H k (t) respectively. 


•j CTJLV 

G k (t) = T 1 T J 

VTn J 


c+iw c k 


ES«. 


2b J E,(sr 

c — ICO 1 V ' 


- oo < t < qo, a < c < b 


(3.45) 


-« c+ioo i~, / \ 

H k (t) = — f 5lWe a ds, 


-oo<t<oo, a < c < b 


(3.46) 


the integral converging absolutely. 


Thus we may consider the simultaneous dual convolution transforms 
(3.40) and (3.41) of secion 3.2.1, where f k (x) and g k (x) are given and <|)h(x) are 
to be determined. 


3.2.3 The Reduction of (3.40) and (3.41) to Transform with A Common 
Kernel 


As usual, we reduce two given transforms with a common kernel. Then 
problem is reduced to inverting the sole convolution transform. To attain this 
object we suppose further that 


j i) £ _i i_ 

j=lv a fj a 3,jJ j=I V a 2,j a 4,j 


converge and that 
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s 


: ^1 -83 +2 


j=iV a hj 


O 

a 3 ,jy 


-8 2 —84 +y 


1 


1 


j=lV a 2 ,j d 4 ,]) 


Now as in proceeding section, if 


00 

f 

E 4 (0 + ix) 

J 

-00 

E| (0 + ix) 


dx < qo for every 0 in a < 0 < b , 


1 


Effcr + ix) 


Ej (0 + ix) 


dx <qo for every 0 in a < 0 < b , 


li m E 4 (0 + ix) 
H 00 e| (0 + ix) 


= 0 uniformaly in a < 0 < b 




lim E| (o + ix) 

— — ; = 0 uniformaly in a < 0 < b , 

T ->oo Ej(0+ix) 


Then 


1 TE 4 (s) 


GUt) = — f st d 

tW W J„E|(s) e 


-oo<t<oo, a<c<b 


(3.47) 


1 CTloo / \ 

H k(t) = — f -"j-r e st dt, - 00 < t < co a < c < b 

2ra cio E l( S ) 


(3.48) 


are well defined and 


I 


Gk(t>- a dt = Ml, 

m) 


a <0 <b, 


(3.49) 



71 







. iJili 












a<cr<b, 


(3.50) 


the integrals converge absolutely. 


We know that G k and H k (t) are frequency functions and that 


Ol(t) 


V 

_q k (t)exp (c4 k t)+ o(exp (a? + s| 


; ( t )f p k (t)exp (aj k t)+ o(exp (a* - s| 


t> O 
t— > -00 


t — ^ +00 

t<n 


(3.51) 


(3.52) 


for s > 0, k - 1, 2, ...., n, where p k (t) are real polynomials and 
a* =max ay, a? min a k j(j = 1, 2, ). 


Since (a, b) is include in Ii n I 2 n I 3 n I 4 , it follows that the bilateral 
Laplace transforms of G k (t) and G k (t) and H k (t) and H k (t) have common 
region of absolute convergence and hence by the product theorem. 


E 4 (s) 


ry = |e^G k .Gj(t)dt= |e _st H k .H k (t)dt 

\ / _on —on 


(3.53) 


both integrals converging absolutely. 


Using the uniqueness theorem of bilateral Laplace transforms, we 


obtain 


G k *G k (x) = H k *H k , k = 1,2, ,n. 


(3.54) 
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for all x. Hence we take the convolution of both sides of (3.40) and G* k (x) 
and of (3.41) and (x) formaly and we have, 


G k f k( x )= J G k(x-u)f k (u)du 


(3.55) 


|Gk(x-u)f k (u)du 


^ w H 

J G* k (x-uJjgkCu-tEa^H (t) dt 


7 n 

= J G k *G* k (x - t)£ (t)dt 


x> X + Q, k= 1, 2, , n. 




fe,! | 


n r 1 n 00 

Z c hkl H k*gk( u )]=Z c bk J H k(x-u^k(u)du 


(3.56) 




S c hk J H k(x-^k(u)du 


II * ** w II 

= Z c hk | H k ( x - u) J H k (u - t)]T b^h (t)dt 
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H k(x-tE a Rk <l>h(t)dt 


where Chk are the elements of the matrix [ahk][bhk] 


Therefore, if we set 


The dual convolution transforms (3.40) and (3.41) can be reduced to be 


transform with common kernel Gk* Gk 


3.2.4 The Solution of (3.4) and (3.41) 


Now if, <j>h(t) is continuous, -oo < t < oo and when (a,b)s 


then on using Widders theorem [129] we can invert equation (3.57) to get 
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- co < t < oo, a < c < b 


c+ioo „w \ 

G k (t) = — f 52^ 


E°k* f k( x )] 


n P / n \ L VJ k* 1 k V x /J 

( t>h(t)=Zd hk - 1 ^ « r * M 1 

h=l ^4(0) 2 c hk [Hj 1 *gj c (x)| 


x > X + Q 
x < A. + 0 




where D stands for differentiation and the operator Ei(D)/E 4 (D) must be 
interpreted as usual meaning, and dhk are the elements of matrix [ahk]" 1 . 


3.2.5 Applications 


To avoid unnecessary confusion we consider 


J HH£ ; 1)2 a^fh (u)du = g x (x), 0 < x 


(3.59) 


.a .. n 


jHH;/; 1 )2b hk f h (u)du = h x (x), 


(3.60) 


k = 1,2, , n. 


In this case, the Mellin transform of H(x)“^ a a ;l)with variables s(= 0 + it) 

coverges if 20 a < (f3 k - a) - 1 and its generating function is T(a + sa)/r(|3 k - 
sa) and for the H-function of (3.60) same results with a replaced by X k and (3 k 
replaced by p hold. 


After an exponential change of variables (3.59) and (3.60) become 
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Hi 


,a ’ 1 ]Z a hk f h( et Vdt = g k (e x ) x> 0 


(3.61) 


k a U 


-M X ’ a . 1 


Z b hk f h( et )e t dt = h k (e x ) x<0 (3.62) 


k= 1,2, 


which are of the form (3.40) and (3.41) with <j>h(t) = f h (e l ) e l 


From these facts, it is easily seen that entire function, defined in 
section 3.2.2 are 


El( S ): 


r(a + sa) 


> E i( s ) = -tr 


\p -sa) 


E 3(»)=zrr ] — v E k (s)=-4 


m K +sa 


r(ji - sa) 


and that 


i; ’ i 


? - 




a + j k P k +j k ^ k +j 


3 J a ’ a 4,j 


k = 1,2,3, ,n and j = 0,1 ,2,3,. 


Here it is assumed that a > 0, J3 k > 0, A. k > a + 1, p > |3 K +1 and from 
familiar infinite expansions of Gamma function, it is clear that Q defined in 
section 3.2.2 is zero. 




After simple calculation, we have 
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G k($ 




t<0 
t > 0 


(3.63) 




o 


t < 0 


,>0 


(3.64) 


where b = 1/a, k = 1, 2, 3, ,n. 


In order to reduce given simultaneous dual equations (3.59) and (3.60) 
to which have a common kernel, we can use two operators of fractional 
integration defined by x and R, given by Fox [41] which after an exponential 
change of variables can be written as 


TW' 


(e~ x )= - p - b - ~J (l - e b ( x - v )^ k - 1 e b p k ( x - v )w(e~ v )lv, 0<x<l (3.65) 


Rw 


^^ = F^^ 1 " eblx " v, ^' a ‘' eb “ <x ' v)w ^ v ^ 


v, x < 1 (3.66) 


Thus from (3.63) and (3.64), it is clear that 


xw(e x )= Jg£(x- v)w(e V )dv = Gjj.*e(e x ) 


(3.67) 


Therefore from (3.67), we have 






.'.i! 


m h , - 




■P . , ^ 


BUMf 
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Gk*g k (e" x )=xg k (e“ x ) 


(3 68 ) 


Similarly from (3.64) and (3.66), we get 


E c hk( H k*h k (e X ))=S c hkR h k(e" X )- 


Hence the solution of (3.61) and (3.62) is given by 


“ “ E 40%„j4H 


x>^+Q 


x<\+Q 


k = 1, 2, n. 


where dhk are the elements of [ah k ]' 


Similarly, we can solve the equations (3.59) and (3.60) of order n by 
operating x and R respectively. 




__n., - 


CHAPTER -4 


In this chapter we have considered the solution of certain triple integral 
equations with Legendre function as kernels by using the properties of 
generalized Legendre functions and inversion theorem for the generalized 
Mehler-Fock transform. This chapter contains brief introduction in section 
4.1, some results in section 4.2, triple integral equations of the first and 
second kind in sections 4.3 and 4.4 respectively, some extensions of triple 
integral equations in section 4.5 and the more general triple integral equations 
in section 4.6 have been given. 


4.1 INTRODUCTION 


Triple integral equations arise in many problems of methematical 
physics [104], Triple integral equations with Legendre functions as kernels 
were first considered by Srivastava [108], He used a method similar to that of 
Cooke [13] to solve the each set of triple integral equations which has been 
reduced to a Fredholm integral equation of the second kind best treated by 
numerical method. Here we have considered more general equations than 
considered by Srivastava [108] or elsewhere. By using the properties of 
generalized Legendre functions and inversion theorem for the generalized 
Mohler-Fock transform, we have been able to reduce the two sets of triple 
integral equations to Fredholm integral equation of the second kind. Further 
extensions to include more general cases have also been discussed. 


i : 



!■ f. 





1^ i 
1 1 ! 
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4.2 SOME RESULTS REQUIRED IN SEQUEL 

It is convenient to list here some results for ready reference. The 
generalized Legendre function P R (z) is defined by [35], 


p£G0= 


r(i-n) 


^Z + l^ 

VZ-ly 


P f 
2 F| 


-v,v + l;l-n;i-izll-z<2 (4.2.1) 
2 2 ) 


The following two integral representations for P^(z) are basic tools for 
our present investigations: 


H-1T 


( 71 ) 

-1/2 


(i Y1 

- 


r 

— u 

v 2 ) 



U JJ 


sinh^a 

J cos (txXcos ha - cos ht) -fi ~2 dt 


(4.2.2) 


where Re(p) < and 


P _i +k ( cosha ) = ( 2% ) ] 


1/2 


f i 


i > 

r|i + p 
1 j 


r — u + ix 

. U j 


p-ix 


-1 


sin h |1 a cosec h(rcx) 

J sin (txXcos ht - coshaX _ 2 dt 


(4.2.3) 


- 1 a 


where - — <Re(p)< — 
2 2 


Let F c (x) and F s (x) denote the Fourier cosine and sine transforms of 
f(x). Then (4.2.2.) can be written as 
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U(a - 1) (cos ha - cos ht ) ^ 2 


= r 


a ^ 


sin h'^aP^, . (cos ha) (4. 2. 4) 


J 


T+« 


where U(x) is the Heavisided’ unit function, so that by the inversion 
formula for the Fourier cosine e transform, we have 


U(a - tXcosha - cos ht) ^ 2 =(2/tc) 1/2 fI 


f\ V 1 


-\x 


sinh ^a 


J 


(4.2.5) 


Rv hk (cos ha)cos (tx) dx 


Also (4.2.3) can be expressed as 


-u/1 ^ 


= tTT 


l 9 +ti 
V2 j 


1 


U(a - tXcosht - cos ha X 2 

sin h^aP 11 , . (cos ha) 

-4+rc v ' 


r — Li+ix 

U H j 


' f\ . ' 


(4.2.6) 


— U - IX 
v2 * j 


which on inversion gives 


U(t-aXcasht-cosha) >i 2 = 2 1/2 tt 3/2 F — + p ]sin hRa 

\1 ) 


Wf-H 


1 


A 


(4.2.7) 


p — ix 


sin (tcx)P^j (cos ha) sin (tx)dt 


J 


If (p(t) is monotonic strictly increasing and differentiable for a<t<b, 
and (p'(t) * 0 in this interval, then the solutions of the equations 


1 


f(x)d> 


[<p(t)-<p(x)f 


M 


actcb, 0<a<l 


(4.2.8) 



* ■ 




[ 1 

■ s 





; i 


i : 


I ' , 

JiLl - „ 

Mxwr. 


^ j 



: 1 


vml 

pi 




-Vxj . V 

1 .'^ 

,u j* mm 
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and 




I 


f(x)dx 


J t kx)-cp(t)r 


: s(tX 


a < t < b, 0<a<l 


(4.2.9) 


are given by. 


f(x)= 7i 1 sin (rax)— f 


<p'(t)g(t)dt 




1-a 


(4.2.10) 


\ -1 • ( \ d f cp'(t)g(t)dt 

f (x) = —it sin (rax) — — r - v /ov 7 


dx J x [ ( p(t)-(p(x)] 


1-a 


(4.2.11) 


respectively. Therefore in view of (3.2.8) and (3.2.10), the integral 
representation (4.2.2) gives 


cos(tx) = (2te) 1/2 F 


(l } 


-p 


cos 


VZ J 




t sinh 1 ^aP^j . (cos ha) 


-t +K 


dt o (cosht-cosha)2 11 


(4.2.12) 


Also, in view of (4.2.9) and (4.2.11), from (4.2.3), we get 


cos 


(tx) = 2 1/2 7C 3/2 rf— + |i) COS (p.7t)rf— - p +ix 


+ sin 


h(rcx)J- 


ooSinh^aP^, . (cos ha)da 


— u - IX 
V2 * j 


(4.2.13) 


~2 +1T 


(cosha-cosht) 2 +>1 


Another result that we shall need is the following inversion formula for 
the generalized Mehler-Fock transform due to Braaksma and Meulehbeld 
[Theorem 7] for conditions of velidity. If 
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H 




4 


! ; ■ a 




: ■ fH 




! * h,‘j 



(4.2.14) 




(p (cosha) = | P^i +ix (cos ha)f (x)dx 


then under certain conditions 




1 


U. + 1X 

V2 * j 


1 


f(x)=tt 1 tsinh(7tx)rj 

00 

J (cos ha)(p (cos ha) sin hada 


\ 


— LX — IT 
v2 * j 


(4.2.15) 


4.3 TRIPLE INTEGRAL EQUATIONS OF THE FIRST KIND 

We with to solve here the triple integral equations of the first kind of 
the form 


J(p(x)rf|-p+ixV • 


i 

--P-1X 


•'s 


12 


J 


P^, . (cosha)dx = 0, 


--HT 


sinh(ttx) 
0<a<aj 


(4.3.1) 




> ■! 


J (p(x)P M ’ 1 (cos ha)dx = f (a), aj < a < a 2 


(4.3.2) 


H 


1 

x U | — p + ix 


1 


A 


— u-ix 

V2 j 


sinh (ttx) 


P ^i l ix (cosha)dx = 0, a 2 < a < ao 


(4.3.3) 


where f(a) is known and <p(x) is to be determined. 


Now we assume that for ai <a <a 2 . 
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(4.3.4) 


00 Tl A f'-i \ 

| cp(xX| — - (i + ix r — [i-ix x sinh (rex) 


ptl i ^.(coshot^x = x(a) 


where %(a) is yet specified. On using the equations (4.2.14) and (4.2. 15) in 
(4.3.4), we have 


cp(x) = tT 1 J x(a') sinha' P^ i+it (cos ha')da' 


(4.3.5) 


Substituting (p(x) in (4.3.2) and inverting the order of integration, we 


f(a)= j x( a 0 L(a 5 a')sinha'da' 


(4.3.6) 


where 


L(a, a') = — Jp^, ^ (cos ha) P^ i+ix (cos ha')dx 


K o 1 


if] \ 

L(a,a') = — T — [i sin h^a sin h^a' 
tt _ \2 J_ 

min (a, a') 

V dt 

| (cos ha - cosht ) ,x+1/2 (cos ha' - cosht )^ 172 


Substituting the value of L(a,a') is equation (3.3.6), we get 


(4.3.7) 


(4.3.8) 








1 1 : 

ir : 
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f(°0 : 


f 1 ^ 
V 2 J 


% 


|x(a') sin h 1+p, a' sin h^a da' 


(4.3.9) 


min(a,a') 

j 


dt 


(cos ha - cosht) R+1 72 (cos ha ' - cosht)^' 


+ 1/2 


On changing the order of integration, we note that 


a 2 min(a,a') a a' a 2 a a a 2 a 

| d(a') jdt = J da'Jdt + | da'j dt = J dt J da' + J dt J da' 

a x 0 0 a 0 aj t 0 0 


On inverting the order of integration (the integral being understood). 
Hence we have (on replacing t by y in last integral). 


dt 


x(a')sinh 1+!J 'a'da' 
+ 1/2 


a, 


71 


(cosha-cosht) |1+1/2 t (cosha'-cosht)^ 
2 

dy 




- f W- j 


sin h^a q (cos ha - coshy ) M ' +1/2 


(4.3.10) 


a 2 

J 

«i 


%(a')sin h 1+ ^a' da' 
(cos ha' - cosht ) |X+1 12 








T 


If we now write that 


c(l ) J xMsmh'^q'da- 


\ (cosha'-cosht) 11 


.+ 1/2 


(4*3.1 1) 


' Mi 

11 , 

* 

t’-j 


and hence equation (4.2.11) gives 
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x(a')sin h 1+fi a' = -— sin [p + 1/2]— , f s^t ^(OjL — (4.3.12) 

K ^a. «' (cosht-cosha')^ 1 


and by equation (4.3.10), we have 


G(t) - ~ sin [tt(ji + 2 )]^“ J 


sin ha da 


^ a, (eos ht - eos ha) 2 ^ 


44-,)Pf(x)_| dy 


sinh J 0 (cos ha- coshy Y +l 11 


°j * 2 x(“ , ) s i n h 1+p a'da' 
a) (cos ha' - coshy ) y+1/2 


On simplifying the last integral final result is 


(4.3.13) 


G(t)= jG(s)K(s,t)ds = F(t) 


aj 


where K(s,t) and F(t) is given by 


(4.3.14) 


w A _ 1 sin htsinhs sin [71(1 + tt )] a ( 3 (coshaj -coshy)2 R dy 

K(s, t)— — __ 77 ~^ J / , v , , \ (4.3.15) 


71 (cosht-cosha^ 


-p J Q (cos ht - coshy Xcos hs - coshy) 


and 


F(t)= [r(l-p)f sin 


7T 


f O 

P + - 
V ^ J 


d | sin h 1_p 'af(a)da 
^ ttl (cos ht - cosha) 2 _|1 


(4.3.16) 


The equation (4.3.14) is a Fredholm integral equation of the second 
kind. The solution of the equation (4.3.1) to (4.3.3) follows from the 
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equations (4.3.5), (4.3.12) and (4.3.14). 


4.4 TRIPLE INTEGRAL EQUATIONS OF THE SECOND KIND 

Now we consider triple integral equations of the second kind of the 

form 


J xt|>(x)P^j (cos ha)dx = 0, 0 < a < ct] 


(4.4.1) 


j(p(x)rf^--[i+ix 


\ 

f 1 . ^ 

r 

--p-ix 

J 

v2 J 


sinh(7tx) 


. (cosha)dx = f(x), aj < a < ct 2 


Hit 


(4.4.2) 


J x(p(x)P^j ^ (cos ha)dx = 0, 


a 2 < a <qo 


(4.4.3) 


where f(a) is known and cp(x) is to be determined. Let us assume that for 
cti < a < a 2 . 


J xcp(x)P^ 1 (cos ha)dx = g(a) 


(4.4.4) 


where g(a) is yet to unspecified function. On using the equations (4.2.14) and 


(4.4.5) 


(p( x )=-r 


it 


1 


\ 


Ll + lX 

V2 . 


if-- H -ill sinh (rex) 


a 2 

J ^-i+h ( cos ^ a ) g(a') sinha'da' 


a i 


Substituting the value of <p(x) in (4.4.2) and inverting the order of 
integration, we have 


a 2 

f(a)= J g(a ') L ! (a , a ')sinh a 'da ' 


where 


1 00 

L(m,h')=— f 

k i 


J 

r — u+ix 
. U j 


\ 


l2 


u 


[X-1X 


J 


sinh 2 (7tx) 


P^ 1 , . (cosha)P^, . (cosha')dx 


2 +1X 




Now using the equations (4.2.3) and (4.2.6), we have 


L(a,a')- 



- 

d Y 

r 2 

% 

r 





v 2 J. 

. 


7 J 


dt 


sin h^a sin h*V max ( a>a ') (cosht - cosha')^ 1 " 
1 


(cosht -cosha)^ ^ 

On using the value of Li(a,a'), we get 


(4.4.6) 


(4.4.7) 


(4.4.8) 


88 



oc 2 , 

f(a) = 7t f — 
Z sin 


sin ha' 


KH 


“i-2 


f 


h^a 


sin h^a' 
dt 


max(a,a')( c °sht-cosha)2 11 (cosht - cosha)2 H 


(4.4.9) 


On changing the order of integration, we note that 


oc 2 co a oo «x 2 oo a 2 t co oc 2 

J da' Jdt = J da' j dt + J da' J dt = Jdt J da' + Jdt J da' 

ccj max (a, a') aj a a a' a 2 otj a 2 04 


On inverting the order of integration. 


Hence we have (on replacing t by y) in last integral. 


dt rg(a')sin h 1 ^a'da' 


- — J 

(cosht - cosha)2 ^(cosht-cosha')* ^ 


-,-2 


-sinh^af(x)~ J 


(coshy - cos ha)* ^ 


°f g(a')sin h 1 **a' da' 
(coshy - cosha')* - * 1 


If we now write that 


J g(ct')sin h 1 ^a' da' 
a , (cosht - cos ha')*"^ 


and hence from (Abel integral equation) (4.2. 1 1), we get 


(4.4.10) 


(4.4.11) 


gMsinh-V^inU I-,]]! f sinllt °M Jt 

71 ’ * 


2 j I dt 


a r 


1 (cosha'-cosht) 2 ' 


Now, by equation (4.4.10), we have 


1 . r (x \ld°r sinhada 

G (t) = _ j — 

71 dt t (cos ha - cos ht)^ ^ 


[fe+^F 


00 J 

sinh 11 a f(x)- f — . 

a2 (coshy - cos ha ) 1 ^ 


°r g(a')sin h 1 ^a'da' 

^ (coshy - cos ha' ') 1/2_M - 

On simplifying the last term, the final result is 

«2 

G(t) + | G(s)Kj (s, t)ds = F(t) 

a, 

where Ki(s,t) and F(t) is given by 

. , sin ht sinhs sin [tc(^ - |i)](cos ha 2 - cosht) 2 ' 


K,(s.t) = 


7t(cosha 2 -coshs) 2 ‘ 


(coshy - coshtXcoshy - cos hs) 


)J 




It 2 

dt !( 


(4.4.12) 


(4.4.13) 


(4.4.14) 


(4.4.15) 


(4.4.16) 
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The equation (4.4.14) is a Fredholm integral equation. Therefore, the 
solution of the equation (4.4.1) to (4.4.3) follows from the equations (4.4.5), 


(4.4.12) and (4.4.14). 


4.5 EXTENSIONS 


If the right hand sides of equations (4.3.1) and (4.3.3) are hj(a) and 
h 2 (a) respectively instead of zero, we may proceed as below. 


First we determine a function A(x) which is such that 


jAWrfi -p + ix) rf — -p + ixlx sinh (rex) 


P^j . (cosha)dx = h] (a), 0 < a < aj 


(4.5.1) 


| A(x)P^, + .^ (cosha)dx = (a), a x < a < a 2 


(4.5.2) 


00 /• | \ N 

JA(x)F — p + ix T — -p-ix xsinh(7rx) 


p V ■ (cos ha)dt = h 2 (a), a 2 <a <oo 


(4.5.3) 



where 1(a) is any convenient function of a. A(x) can be determined by 
inversion from (4.2.15) and so we can find a function m(a) which is such that 


J A ( x fc +k ( c osha)dx = m(a), for < a < a 2 


(4.5.4) 






Now write B(x) = <p(x)~ A(x)and subtract equations (4.5.1), (4.5.4) and 
(4.5.3) from equations (4.3.1), (4.3.2) and (4.3.3) in their now form and we 


have 


Ml 


\ 

ft + ir 
2 ) 


jj. — ir I x sin: 


inh (tit) 


P^i + j T ( cos ha)dx = 0 0<a< ai 


(4.5.5) 


| B(x)P^ 1+ .^ (cosha)dx = f (a)- m(a), aj < a < a 2 (4.5.6) 

o 2 


Ml 


g + ix 


■(J.-1X X SI: 


influx) 


( cos ha)dt = 0, 0 < a < qo 


(4.5.7) 


which are of the form (4.3.1), (4.3.2) and (4.3.3) can thus be found and so 
(p(x) is found. 

4.6 THE MORE GENERAL TRIPLE INTEGRAL EQUATIONS 

The equations considered in section 4.3 can be further generalized if 
there is an extra factor [1 + H(x)] in equation (4.3.2) in its now form, where 
H(x) is a known function. Thus case may be treated in the same case as done 
by Cooke [21] but it loads to multiple integrals, which are not specified. Then 
the final result [5.3.14] will be a Fredholm integral equation of the second 
kind with an extra term for the new form of the second kind with an extra 
terms for the new form of the equation. The extra term will be 
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J G(s)ds J H(x)l(y, x)l(y, x)d 


( 4 . 6 . 1 ) 


where 


( 4 . 6 . 2 ) 


CHAPTER -5 



TRIFLE SERIES EQUATIONS INVOLVING ASSOCIATED I 

LEGENDRE FUNCTION ; II 


5. INTRODUCTION 

Triple series equations are useful in finding the solutions of three part mixed 
boundary value problems of electronics, elasticity and other field of mathematical 
physics [104] Collins [1 1] solved the triple series equations for the first time and since 
then many authors like Lowndes [55], Dwivedi and Gupta [27], Parihar [78], Melrose 
and Tweed [62] obtained the solutions of the triple series equations involving types of 
special functions. In the next section we find solution of some triple series equations 
involving associated Legendre function. 


5.1 THE EQUATIONS 

Here we are concerned with triple series equations of the form: 
(i) Triple Series Equations of the First Kind 

CO 

S(2n + 2m + lK«T;"(Cose) = f 1 (et O£0<a 

n=0 


CO 

Z AnT»«(Cos0)= f 2 (ei a < 0 < p 

11=0 
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3 < 0 < Tt 


(5.5.3) 


Z (2" + 2m + lK T‘“ (Cos0) ={,(&), 


n=G 


(ii) Triple Series Equations of the Second Kind 


Z B nT^„(Cose)= gl (9)i 

n=0 


0 < 0 < a 


(5.5.4) 


Z (2n + 2m + l)B n T m ” (Cos0) = g 2 (G), 

n=0 


a < 0 < p 


(5.5.5) 


Z B n Tm?n(Cose) = g 3 (e) 


n=0 


3 < G < 71 


(5.5.6) 


i* 


In above equations fi(G) and gi(0), (i = 1,2,3) are the prescribed functions of the 
variable 0 and the equations (5.5.1) to (5.5.6) are to be solved for the unknown 
coefficients A n and B n . T m “ (Cos 0) is the associated legendre function of degree n+m 
and order -m of the first kind, 

5.2 PRELIMINARY RESULTS 


■ -..fm 


(i) Inversion Theorem for Associated Legendre Polynomials. If We Assume That 
The Expansion 


f (0) = Z (2" + 2m + iJC.T"" (Cos0) 



is valid for 0 < 0 < tc and that it can be integrated term by term, the coefficient C„ are 
given by 


c„ = 2(-)”f f(xfc„(Cosx)Sinxdx 


(5.5.2) 


The Series 


S m (6, x)=2 T (Cos0)T“ +n (Co 

n=0 


(5.5.3) 


/_\m oo 2 

osr)Cosmv|/d\j/ 


(5.5.4) 


2(~) m j-min(s„s 2 ) 

7t( Sl S 2 ) mJo I 


s 2m ds 


L 2 2 Y 2 „ 

t?l _s A s 2 ~ s 


2_,2 W2 


(5.5.5) 


Where s, = 2Sin— Cos— x, 
2 2 


So =2Sin-xCos-0 
2 2 


And both Si and s 2 are positive for all 0 and x since 0 < 0 , x < it. 


and S = 2Cos— Cos— tan— 

2 2 2 









(iii) Finally, We Require An Integral Representation of the Associated Legendre 


Polynomial T m+n (Cos 8), to Obtain Which We Make Use of the Following 


Result: 


If, f(9)=j! 


(Cosu - Cos0) /2 


“ d f '( e >=£rr 1 


(Cos0 - Cosu) i 


a < 0 < (3 


a < 0 < p 


(5.5.6) 


(5.5.7) 


then - de , 

’ rduJt ‘(Cos8-Cosu)>i 


< 16 . 

Itdu “(Cosu-Cose)/^ 


a < u < P 


a < u < P 


(5.5.8) 


(5.5.9) 


■B 


THE SOLUTION 


5.3 EQUATIONS OF THE FIRST KIND 


In order to solve the triple series equations of the first kind, we set 


Z(2n + 2m + l)A n T-“(Cos0) = h(0), a < 0 < p (5.3.1) 


Applying the inversion theorem (5.2.2) in equations (5.1.1), (5.1.3) and (5.3.1), 


we get 
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A " = [f f l ( X ) + £ h ( X ) + £<3 ( X )] T »« (Cosx)Sinxdx (5.3.2) 


Substituting the expression (5.3.2) for An in equation (5.1 .2), we get 


£T-J„(Cose)..L(-)” 

n=0 1 




: T™ +n (Cosx)Sinxdx = f 2 (0), 


a< 0 <p 


(5.3.3) 


Interchanging the order of summation and integration, we get 


2 




Z T m+n (Cos0)r^ +n (Cosx)x Sinxdx = f 2 (0), 

n=0 


a < 0 <P 


(5.3.4) 


Using the expression given by (5.2.3) in (5.3.4), we obtain 


(-)m h(x)S m (0, x)Sinxdx =F(0), 


a < 0 < P 


(5.3.5) 


where 


F (0) = f 2 (0)-| (-) m fj (x)S m (0, x)Sinxdx + £f 3 (x)S m (0, x)Sinxdx 


(5.3.6) 


Now using the summation result in terms of integral (5.2.5) in the equation 
(5.3.5), we get 
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9 r 

j y x) r»<^) — g * . Sinxdx 

*(^ 2 ) |, 2 _ s 2 )( S2 2 _ s 2|/2 


= F(eX 


a < e < p 


(5.3.7) 


After breaking into parts, and noting that si > s 2 when 0 < x, and si > s 2 when 


0 > x, we have 


f h(x)sin x f ™ dx 

* 1 J “ 


s 2m ds 


f 1 1 1 r 

. „ 2 cos - 0. cos -x. tan— u 

— — - f h(x)sin'- m xdx.— — < — du 

71 “ 2 0 [(cosu - cosOXcosu - cosx)] /2 




+ j^h(x)sin 1 m xdx. 


fill Y 
2cos-0.cos-x. tan-u 

i£A_? ? 

^ [(cos u - COS ©Xcosu -cos x)]/2 


= F(0), a < 0 < P 


(5.3.9) 


j>) 


sin 1_m x f l 


s 2m ds 


" s 2 I S 2 2 -S 2 )} 2 


Tdx = F(0), a < 0 < P (5.3.8) 


Making change of variables by putting the values of s, si and s 2 , we get 


, 1 , 5 I 









Now inverting the order of integration, we get 




Interchanging the order of summation and integration, we get 


2 


f)"[f k ,(x) + j!g 2 (x) + j^k 2 (x)]|;T m ”(Cos6) 


n=0 


Tm+n (Cosx)x Sinxdx = \ 


ate). 

gate) 


O<0 <a 
P < 0 < 7t 


Using the result (5.2.3) in equations (5.4.6) and (5.4.7), we have 


( X ) + Jp k 2 (*)}s m (e,x)Sinxdx = | 


(5.4.6) 

(5.4.7) 


M(0), O<0<a (5.4.8) 

N(0), P < 0 < 7t (5.4.9) 


where 


M ( e ) = 81 (®) ■ ~ ~ (~) m £ g 2 (x)S m (0, x)Sinxdx 


and N(0) = g 3 (0) - - (-) m f g 2 (x)S m (0, x)Sinxdx 


(5.4.10) 

(5.4.11) 


Now using the summation result (5.4.5) in terms of integral in equation (5.4.8), 
we obtain 


1 ( \2m 


2 . — (-)' 
2 V 1 


tt( Sl S 2 f 


r^iwr 


(s,s 2 ) S 2m ds 


B , 2 -s 2 )s 2 2 -s 2 f 2 


-.Sinxdx 


+ 


/ \ r mmi 

p k 'Wi, 

= m(o). 


(s,s 2 ) S 2m ds 


s> 2 - s 2 (s 2 2 -s 2 |' 


-Sinxdx 


0<9 <a 


(5.4.12) 
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Sin ~ m 0 


J^k^x^in 1 m xdxj^ 


S 2m ds 

[s, 2 -s 2 )(s 2 2 -.s 2 )} ! ' 


J^k^xJSin 1 ”xdxj^ 


S 2m d$ 


is, 2 -s 2 X s 2 2 - s 2 )} V; 


■ j\ 2 (x)Sin 1 m xdx[ S ' 


S 2m ds 


is , 2 -s 2 Xs 2 2 - S 2 )]' 2 


= M(e) 


0 < 0 < a 


(5.4.13) 


Now changing the variables by putting the values of Si, S 2 and S in equation 


(5.4.14), we get 


( r Q r x u'] 

2Cos— Cos-tan — 

2=22 fk 1 (x)sin 1 -”xdxij' s i 2 - 2 ^ Hu. 

% ] 2 * [(Cosu - CosOXCosu - Cosx)]/2 


( Q x uV 
2Cos-Cos — tan — 

+ r k j (x)Sin 1_m xdx - — — du 

2 [(Cosu - CosGXCosu - Cosx)]/ 2 






0„ x uY 

■ f*k 2 (x)si„ ‘-“xdx Y— i 22 Y^- du 


^ 0 [(Cosu - CosOXCosu - Cosx )]/ "2 


M(et 


0 < 9 < a 


(5.4.14) 


f > ' : . Y 


• P • 


. . ... - 


Inverting the order of integration, in the above equation, we have 


2~ m Sin ~ m 9/2 Cos~ m 6/2 ,2 2m Cos 2m 6/2 T re tan 2m u/2 

271 * (Cosu - Cos Qf 2 ^ 


) 2 m k t (x)Cot m x/2 .Sinx re tan 2m u/2 
(Cosu - Cosx) 1 2 0 (Cosu - Cos0) 1/2 

* 2~ m ki(x)Cot m x/2 .Sinx + r0 tan 2m u/2 
(Cosu - Cosx) 1 ' 2 '° (Cosu- Cos6) 1/2 * 


r82 kj^xJCot x/2. Sinx r9 tan" 111 u/2 , 

■ Ju //Z ~n V/2 dX + JO ZTZwJ 411 


rg 2 m k 2 (x )Cot m x ' 2 .Sinx 
P (Cosu - Cosx) 1 2 


M(e), 


0 < 0 < a 


2 2m tan~ m 0 2Cot 2m 6 2 r 8 tan 2m u/2du 
2tt _ 0 (Cosu - CosQ) 1 ' 2 

f<x k{(x)Cot m x 2. Sinx ^ + r8 tan 2lB u/2 
(Cosu -Cosx) 1 2 ' 0 (Cosu-Cos6) 1/2 U 


tk 2 (x)Cot m x 2 .Sinx 
(Cosu - Cosx) 1 2 


M(e), 


0 < 0 < a 


Now above equation can be written as 


1 rQ tan 2m u 2 ^ f® k t (x)Cot m x/2 .Sinx ^ 

2tc ■*» (Cosu - CosG) 12 U J u (Cosu - Cosx) 1/2 


+ pk,(x)Cofx 2.Sinx dx = ^ I ftM(e) , 0 £ 0 < , 
P (Cosu - Cosx) 1 ~ J 2 


Equation (3.8.16) is now reduced to the following form 


Mr- 


(5.4.15) 


(5.4.16) 


I 

I 




ifc 


■IS Bt 1 . : 
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(5.4.17) 


£■ 


tan*" u 2 


(Cosu - Cos©) 1 


J/2 



k 2 (x)Cot m x/2 .Sinx 
(Cosu - Cosx) 1 2 


■dx 


= tan m — 0.M(0), 
2 


0 < 0 < a 


where 



kj(x)Cot m x/2. Sinx , 

— (jx 

(Cosu - Cosx) 1 2 


(5.4.18) 


Solution of the equation (5.4.17) can be given, with the help of the results (5.2.6) 
and (5.2.8), as 


j/ \ J_ j- rck 2 (x)Cot m l 2x,Sinx _ Cot 2m u/2 
U + (Cosu -Cosx) 1 2 X_ ic 


d p M(0)tan 2m 1 29.SinQ 
du 0 (Cos0 - Cosu) 1 2 


0 < u < a 


(5.4.19) 


Now equation (5.4.19) becomes 


\ % 1 fwk 2 (x)Cot m x/2 Sinx . 

M,(u) = J(u)+ — L-4 1 - z : — -dx, 

P (Cosu-Cosx) 1 


0 < u < a 


(5.4.20) 


where 


M| (u) = 


Cot 2m u/2 d f uM(e)tan m l/29.Sin6 


ic du*> (Cos0 - Cosu) 12 


d0 


(5.4.21) 





Again starting from equation (5.4.9), we have 


(-) m f k i ( x )+ Jk 2 (x) S m (6, x)Sinxdx = N(0), p < 0 < k 


Using summation result in terms of integral by equation (5.2.5), we get 


1 / \2m 


2 . — (-)' 
2 V ' 




(SiS 2 f 




■mm! 
'0 


in(Si,S 2 ) • S 2m ds 


(s, 2 -s 2 )s 2 2 -s 2 f- 


-.Sinxdx 


+ 


fMr 


m(s,.S 2 ) 


S 2m ds 


(s, 2 -s 2 )s 2 ~s 2 )]" 


-Sinxdx 


N(e) p<9<7C 


Sin~ m e 


71 


k ! (x)Sin 1_m xdx£ 2 


S 2m ds 


^ 2 -s 2 ^ 2 -s 2 )F 2 


+ £k 2 (x)Sin'-”xdxf 


S“ m ds 


'(s, 2 -s 2 | 2 2 -s 2 )f 2 . 


= N(e), P < 0 < k (5.4.22) 


Now making change of variable in equation (5.4.22) 


Sin _m 0 


71 


£k,(x)Sin I - m xdx.-j* 


1 rx (2Cos 0/2 Cos x/2 tan u/2)' 


,2m 


2 J ° [(Cosu - CosOXCosu - Cosx)] - 2 


du. 


fk 2 (x)Sin 1 “ m xdx.-t [ x (2Cos8.2.Cos^tg j ^" 

^ 2 0 [(Cosu - CosOXCosu - Cosx)]/2 


+ Ck.txJSin'-xdx if (2C°se 2.COSX 2 2 tm m2 ) ^_ du 
2 [(Cosu - Cos0XCosu - Cosx)]/2 

= N(0> P < 0 < % 


(5.4.23) 






no 



Inverting the order of integration in equation (3.4.23), we get 


2~ m Sin “ m 0/2 Cos -m 6/2 .2 2m Cos 2m 0/2 f r a tan 2m u/2 


(Cosu - Cos©) 1 2 


2~ m k 1 (x)Cot m x / 2.Sinx 
(Cosu - Cosx) 1 2 


, rP tan u/2 , 

h (Cosu -Cose) 1 '' 2 


•7 t2 m k 2 (x)Cot m x 2,Sinx^ __ + re tan 2m u/2 
P (Cosu - Cosx) 1 2 ■’P (Cosu - Cos0) 1/2 


fit 2 k->(x)Cot x 2.Sinx , ,, 

: — = — dx = N(0), 

u (Cosu - Cosx) “ 


_1_ re tan 2m 12 udu r^ k 2 (x)Cot m l/2x.Sinx^ 
2 ti Jp (Cosu - Cos©) 1 2 “ (Cosu - Cosx) 1 2 


= tan m -0.N(0)-— f tan 12u d u 
2 27t Jo (Cosu-Cos©) 12 

k,(x)Cot m -x.Sinx tan 2m -u 

x f 1 —7 — — rjy-dx f — — -du 

Ju (Cosu - Cosx) 1 2 2k Jo (Cosu - Cos©) 1 2 


(• ^k 2 (x)Cot m 1 2x.Sinx 
P (Cosu - Cosx) 1 2 


P < 0 < K 


Equation (3.4.25) now becomes 


(Cosu - Cos©) 1 2 


^-r-rdu = tan m ■ 


r«-N(e)-r 


j(u)tan 2m — u 

2 

(Cosu - Cos©) 1 ' 


. tan 2m -u k 7 (x)Cot m — x.Sinx 

± (f Z — — d U x r-4— 2-rr-^P < 0 < n 

2 k * (Cosu - Cos©) 1 2 J P (Cosu - Cosx^ 2 





where 


, k 2 (x)Cot m - x.Sinx 


(5.4.27) 


Now putting the value of J(u) from (5.4.20) in equation (5.4.26), we get 


K(u)tan 2m -u 
f° 11 2 

P (Cosu - CosO) 1 ■ 


-rrdu = tan — 

i 12 2 


9.N(e)-r 


tan 2m -u 


(Cosu - Cose ) 1 


Mj(u)- — r 
1V ’ O-n-Jf 


,k 2 (x)Cot m - x.Sinx 


2?t J P (Cosu-Cosx) 1 




j „ tan 2m -u k 2 (x)Cot m -x.Sinx 

— JlT 2 d uf 2 dx,p < 0 < re 

2n * (Cosu - CosG) 1 2 J P (Cosu - Cosx) 2 


(5.4.28) 


„ K(u)tan 2m -u , ,, , u zm . , _ 

j? 7 ^rydu = tan m -0.N(e)-f ! ^ tan — ?™du 

(Cosu -CosG) 1 2 9 20 if '— ^- Ml/2 


ra M| (u)tan 2m 1/2 u 
'° (Cosu - CosG) 1 ' 2 C 


. 2 in 1 

tan - u 


1 fP 2 ~ , r«ko (x)Cot m 1/2 x.Sinx D ^ A ^, 

_ - du | 2— <ix, p<0<7i 

2ft Ja (Cosu - Cos0) 12 J P (Cosu - Cosx) 12 


(5.4.29) 


In view of the results (5.2.6) (5.2.8), equation (5.4.29) gives 




; ! ; ' ■ 
,V ! ,1," '4,* r , 

J " ji 1 | I ' i* 






fa M 1 (a)tan 2 “ia ( 
^ (Cosd - Cosx) 12 da 


(5.4.33) 



Changing the order of integration of the equation (5.4.3), we get 


K(u)=N 1 (u)+M 2 (u)- 


Cot m — u 


— f P tan 2m -ddd 


2 k z L Ja 2 




_d_ ru SinSdS k 2 (x)Cot“-x.Sinx 

du Jp (CosS - Cose) 1 2 (Cose - Cosu) 1 2 ■’P (Cose - Cosx) 1 ' 2 X ’ 


3 < U < 7t 


(5.4.34) 


using the result given by 


Sinede 


du J p (Cose - Cose) 1 2 (Cose - Cosu) 1 2 


Sinu(Cose - CosP) 1 2 
(Cosj) - Cosu) 1 2 (Cose - Cosu) 


(5.4.35) 


In equation (5.4.34), we get 


Cot 2m — uSinu tan 2m — e(Cose - CosfB) 1 2 

K(u)= N,(u)+M 2 (u) — - rrr- f j- — r dv 

V W 1 2k (Cos|3 - Cosu) J “ (Cose -Cosu) 




1 




■ 

m 


m 




l ''' \ 



Now with the help of the results (5.2.7) and (5.2.9), the solution of the equation 
(5.4.27) is given by 


By (3.8.37) we obtain 


K(s)Sinsds 


Putting the value form (5.4.38) in equation (5.4.36), we get 


uSinu 


k(s)Sinsds 


Changing the order of integration of equation (5.4.39), we get 
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Cot 2 -u 

K(u) = N 1 (u)+M 2 (u) - 2 - 


1 k(s).Sinsds 

' (Cosp - Cosu) 1 2 (CosP - Coss) 1 " 




(CosS - CossXCosS - Coss)’ 


tan 2m — SdS 


; (Cos& - CossXCosS - Coss)’ 


P < U < 71 


(5.4.40) 


Equation (5.4.40) can now be rewritten as 


k(u)+ — fk (s)R.(s, u)ds = Nj (u)+ M 2 (u) P < u < n 


(5.4.41) 


where R(s, u) is symmetric kernel. 


R(s,u) 


2m -u.Sins.Sinu 


(CosP - Coss) 1 2 (Cosp - Cosu) 1/2 


(5.4.42) 


s ( s - u )= £ 


tan 2m -SdS 


(Cos& - CossXCosd - Cosu) 


(5.4.43) 


Equation (5.4.41) is a Fredholm integral equation of the second kind which 
determines K(u) and from equation (5.4.37) k 2 (x) can be find. After that, we can 
calculate the value of J(u) from equation (5.4.20) and consequently ki(x) can be found 
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with the help of equation (5.4. 18). Finally, the unknown coefficients B n can be computed 
with the help of the equation (5.3.2) which satisfy the equations (5.1.4), (5.1.5) and 
(5.1.6). 


PARTICULAR CASE 


If we let 3 — 7i, the equation from (5.1.4) to (5.1.6) reduce to the dual series 
equations and the solution obtained here agree with that obtained earlier by Collins [10], 
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CHAPTER -6 

FIVE SERIES EQUATION 

6. INTRODUCTION 

There has been a lot of work on dual, triple and quadruple series 
equations involving different polynomials. Due to the importance of these 
series in finding the solutions of various mixed boundary value problems of 
elasticity, electrostatics and other fields of mathematical physics, a number of 
researchers took interest in finding the series solution as well as developing 
and investigating new classes of series equations. 

There was almost no research work on five series equations until 
Dwivedi and Shukla [32] taken it into consideration. They solved certain five 
series equations involving generalized Bateman K-functions, series of Jacobi 
and Laguerre and the product of V generalized Bateman K-function. In the 
subsequent years Dwivedi and Singh [33], obtained the solution of five series 
equations involving generalized Bateman K-function and Jacobi polynomials 
respectively. 

In the present chapter, we have considered five series equations 
involving series of Jacobi polynomials, which were untouched till date. 

6.1 FIVE SERIES EQUATIONS INVOLVING SERIES OF JACOBI 

POLYNOMIALS 

The solution of five series equations involving series of Jacobi 
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polynomials is obtained by reducing them to Fredholm integral equations of 
the second kind in one independent variable. Dual, triple, quadruple and five 
series equations involving series of Jacobi polynomials can be change into 
the series involving ultraspherical polynomials or Fourier cosine series by 
small amendment to the original one. These latter series equations play an 
important role in solving the mixed boundary value problems, when we 
consider the distribution of stresses in the interior of an infinitely long strips 
containing three Griffith cracks situated on a line perpendicular to the 
boundary lines of the strip. 

Here we are concerned only with five series equations involving series 
of Jacobi polynomials which are extensions of quadruple series considered 
by Dwivedi and Singh. 

6.2 THE EQUATIONS 


We shall solve the following set of five series equations 


co 


z 

11=0 


z 


n=0 



■f.W 

0 < 0 < a 

(6.2.1) 

/ W pM)(cose) = < 

r(n + a + l)r(n+p + 3'2) ' 


b < 0 < c 

(6.2.2) 


f 5 (el 

d < 0 < K 

(6.2.3) 

/ w p^P)(Cos0)=< 

r(n + p + l)f(n + a + 12) V ' 

f 2 (ei 

a < 0 < b 

c<0 <d 

(6.2.4) 

(6.2.5) 


where a, {5 > - V 2 , and f; (0), (i = 1,2,3, 4,5) are prescribed functions and 
equations (6.2. 1) to (6.2.5) are to be solved for unknown co-efficients A„. It is 
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assumed that series (6.2.1) to (6.2.5) are uniformly convergent and f; (0) and 
their derivatives are continuous. 


6.3 PRELIMINARY RESULTS 

In the course of analysis, we require following results: 
(i) The Orthogonality Relation for Jacobi Polynomials 

2a 

J”(Sin 6/2)” (Cos0/2) 2|J p( a ’P)(Cose)pi a ’P)(Cose)Sin0de 
_ r(n + a + l)r(n + J3 + l) 

_ -T^p] 6 mn 

Mn 

Is valid for a > -1, P > -1 
Where 6 mn is Kronecker delta and 

q(cc+P) = i. n t{(a + P + 2n + l)}r(a + p + n + l) 



(6.3.1) 


(ii) The Series 

s( u q \ = y = ql a ' p)r ( n+a+ i 2 ) 

h {r(n + a + 1)} 2 r(n + p + 3/2) 

(Sin u/2f a Pi a ’ P) (Cosu)x p( a ’%os0) (6.3.2) 

(Sin 0 2)~ 2a rmin(u,0) E(y)dy 

it •’ 0 (Cosy - Cosu) 1 2 (Cosy - Cos ©) 12 
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where 


E(t) = (Sin t/2) 2a (Cos t/2)' 2p , t = min (u, 0). 


(iii) We shall use the following two forms of Schlomilch’s Integral 

Equations: 


If f(0) and E(0) are continuous in a < 0 < b, then the solutions of the 
integral equations: 


f(e)= f - — 

a (Cosu - Cos©/ 2 


(6.3.4) 


and f'(0)= J 


b g'(u)du 
3 (Cos0 - Cosu) 1 2 


(u) - — — f u f(6)Sin0d0 
8 U tcdu-k (Cos0 - Cosu) 1 2 


(6.3.5) 


(6.3.6) 


and g'( u ) = -1 A f b _i^)gin9d e 
ttdu •'“(Cosu -Cos©) 1 2 


(6.3.7) 


respectively. 


6.4 THE SOLUTION 


Let us suppose 


g(0), 0 < 0 < a (6.4. l) 

^2 n(a,P)(Cos0) = < h(©^ b <0 <c (6.4.2) 

^r(n+p + l)T(n + a + V2) v ' 

k(0), d < 0 < k (6.4.3) 
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Using orthogonality relation (6.3.1) in equations (6.2.4), (6.2.5) and 
(6.4.1) to (6.4.3) we get 


r(n+q + l.'2) fam 
” r(n + a + 1) q ” 


4 - 


J 0 8 '(u )+ 1 f '2 (u)+ J b C h'(u)x f c f' 4 (u) 

CmSin u/2) ?cx P^ a ’ p )(Cos0)x P^ a ’P)(Cosu)du (6.4.4) 


where 


g'(u) = (Cos u/2) 2p Sinu g(u), 
f 2 ' (u) = (Cos u/2 ) 2p Sinu f 2 (u), etc. 

Substituting the expression for A„ from (6.4.4) in equations (6.2.1) to 
(6.2.3), we obtain 


qjf ,p 4(n + a + 1 2) 


n t?){r(n+a + l)} 2 r(n + p + 3 2)1 

+ jJk’Cu) (Sinu/2) 2a pi a ’^(Cose)P 1 S a ’ p )(C osu)du 


J o a g’(u)+ J a bf, 2 (u)+ J b °h'(u)x £ d f’ 4 (u) 


fi(e) 

0 < 0 < a 

(6.4.5) 

f 3 (e) 

b < 0 < c 

(6.4.6) 


d < 0 < % 

(6.4.7) 


Applying the summation result (6.3.2) and interchanging the order of 
integration and summation, we get 
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where 


[P(e) 

0 < 0 < a 

(6.4.8) 

j^g'(u)+£h'(u)+jjk'(u) S(u,e)du. Q(eX 

b < 0 < c 

(6.4.9) 

We). 

d < 0 < 7i 

(6.4.10) 


p(e)=fi (»)-[£ f'2<u)+ £ f' 4 (u)js(u,e)du 
Q(e) = f 3 (e)- [ £ b f' 2 (u)+ ff 4 (u)l s(u,e)du 


R(e)=f s (e)- | a f 2 (u)+| f 4 (u) s(u,e)du 


Now using the summation result (6.3.3) in terms of integral in equation 
(6.4.8) we obtain 


f a ,, v f mm(u,0) E(y)dy f . bV 2 ” /.x 

Jo 8 u u Jo (Cosy-CosuPfCosy-Cosey 2 ”*1 2 j 


- £h' ( u)du| ) ,nm(u - 6) _ 

(Cosy - Cosu )2 (Cosy - Cosu)! 


BiS 

f ; * ‘i. 


■j;k'(u)duf" <U - 6) 5® T 0 < 8 < a (6.4.11) 

(Cosy - Cosu )2 (Cosy - Cosu )2 


E(y)dy 


fVudu f U - - 

0 0 (Cosy-Cosu) 1 ^(Cosy-Cos©) 1 2 


•jV(u)du| 


E(y)dy 


0 (Cosy - CosupfCosy - CosG) 172 V 2 


n Sin — P(0) 
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rV( U )du f 

Jb w Jo 1 1 

(Cosy - Cosuj2 (Cosy - Cosu)2 


|\'( u ) du ^- i&te . 

(Cosy - Cosu)2 (Cosy - Cosu)2 


0 < 0 < a (6.4.12) 


Changing the order of integration in the above equation, we get 


f 6 E(y)dy r 6 g'( u )d 


f g yujau . r 

■»y / J Jc 


■e E(y)dy 


(Cosy - Cos0)2 (Cosy - Cosu)2 (Cosy - Cos0) 


g f (u)di 


(Cosy - Cosu)2 


Sin- P(0) 


r 9 E (y)dy 


h'(u)du 


(Cosy - Cos0)i (Cosy - Cosu)2 


E(y)dy 


k'(u)du 


(Cosy-Cos0)2 (Cosy-Cosu)2 


0 < 0 < a (6.4.13) 


f E(y»y P S'(u>»u = J Sin «fp(e) 

Jo / J J y / 1 l 2 j w 

(Cosy - Cos 8)2 (Cosy - Cosu)2 v J 


j Q E(y)dy jc h'(u)du 
(Cosy - Cos0)2 (Cosy - Cosu)2 


r® E(y)dy rx k'(u)du 

. j - TJd “ T 

(Cosy - Cos0)2 (Cosy - C0SUJ2 


k'(u)du 


0 < 0 < a (6.4.14) 


"1 
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Using results (6.3.4) and (6.3.6) in eq. (6.4.14), we obtain 


B(y>r 


g'CuJdu 


1H fsioif 

*1-1 r»L V 

7T. flv JO 


P(9)Sined0 


(Cosy - Cosu) 2 ndy ° (Cos6 - Cosy)| 


1 d ry 


l 


SinGdG 


°(Cose-Co S y)j 


1 


,e E(y)dy 


rf- 


h'(u)du 


[ (Cosy - Cos0) 2 (Cosy - Cosu); 




t 

Jo 


E(y)dy k'(u)du 


tI 


(Cosy - Cos0) 2 d (Cosy - Cosu) 2 J 


0 <0 < a (6.4.15) 


Changing the order of integration of above equation, we get 


E(y)j; 


g'(u)du 


f Q\ 2a 

d [ Sin — J p(0)Sin0d0 
~lo 


(Cosy- Cosu ) 2 dy (Cos0 - Cosy) 


K 


f 


h'(u)du d ry 


(Cosy- Cosu ) 2 ^ 


fE(y,dy}; 


Sin0d0 


l 


(Cos0 - Cosy) 2 (Cost - Cos0) 2 


+ 


Jc 


n k'(u)du d ry 


(Cosy - Cosu) 2 dy 


J 0 'E(y)dyf 


Sin0d0 


1 1 

(Cos0 - Cosy ) 2 (Cost - Cos0) 2 J 


0 < 0 < a (6.4.16) 


Using the result 


J, 


Sin0d0 


(Cos0 - Cosy) 2 (Cost - Cos0) 2 


= TC 


(6.4.17) 


h a 
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in equation (6.4.16), we get 


E(y)f 


g'(u)du 


(Cosy - Cosu )'2 
r 21 E(y)k'(u)du 

Jd . 1 

(Cosy — Cosu)! 


r= p i(y)-£ 


E(y)h'(u)du 
(Cosy - Cosu )2 


0<9<a (6.4.18) 


where 




. Sin - 


p(e)sin0de 


(Cos0 - Cosy )2 


(6.4.19) 


Again using the results (6.3.5) and (6.3.7) in equation (6.4.10), we get 
E(y)g'(u) = f a ginyPj^yjdy + _l__d_ p E(y)Sinydy 

7i du •' u / _ vl % du •' u /_ _ vi 

(Cosu - Cosy )2 (Cosu - Cosyja 


+ 


£ 


h'(v)d% 


4 


(Cosy - Cosv): 

k'(v)dv 
(Cosy-Cosv);> 


1 d ra E(y)Sinydy 

i Ju 


1 n du Ju - - 1 


(Cosu - Cosy)o 


0 < u < a (6.4.20) 


Now changing the order of integration and using the result 


d ru Sinydy Sinu (Cosv - Cosa )2 

Jgt 1 1 1 1 

(Cosu - Cosy)'! (Cosy - Cosv)! (Cosa - Cosu)! (Cosv - Cosu)! 


(6.4.21) 
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in equation (6.4.20), becomes 


E(y M u ) = p 2 ( u ) + - £h'(v) A (v, y)dv + jj k'(v)A(v, y)dv 


where 


P 2 ( U) = _i A p_ S nff p i (y)dy 

71 (Cosu - Cosy )2 


0 < u < a (6.4.22) 


(6.4.23) 


A(v,y) 


E(y)Sinu(Cosv - Cosa )2 
(Cosa - Cosu )2 (Cosv - Cosu) 


(6.4.24) 


L 


Again using summation result (6.3.3), in terms of integral in equation 
(6.4.9), we obtain 


Iff. 


h'(u) + rh'(u)ir°' e> 7 — — tit 

6 -I 0 (Cosy - CosuXCosy - Cos0| “ 


^Sin|j Q(0)- £g'(u)du+J\'(u)du 


rm 

X 

Jo 


>min(u,0) 


E(y)dy 


(Cosy - CosuXCosy - CosG/ 2 


b<6<c (6.4.25) 
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j>H"- ^ r 

(Cosy - Cosu )2 (Cosy - Cos0)2 


r 

(Cosy - Cosu )2 (Cosy - Cos0)2 


Sin^l -Q(0)- fV(u)duf U ^)dy 

2 J 7 Jo &v ^ Jo i l 

(Cosy - Cosu )2 (Cosy - Cos0)2 


■ f* k-(u)du f° SfMy 

Jd v ^ JO . 1 1 

(Cosy - Cosuj 2 (Cosy - Cos0)i 


b < 0 < c (6.4.26) 


Changing the order of above equation, we get 


■C - ):K , -j; h,(u)dU , - 4 Sin ffx#) 

(Cosy - Cos0)2 (Cosy - Cosu )2 


jb E(y)dy jc ti(u)du 
(Cosy - Cos0)2 (Cosy - Cosu )2 


h'(u)di 


. j a E (y)dy j a g \ ujgu 

(Cosy - Cos0)2 ’ (Cosy - Cosu )2 


g f (u)du 




> 

I 

| * 


: : 


rQ E(y)dy r* k'(u)du , A 

J®: — b<0<c (6A27) 

(Cosy - Cos0j2 (Cosy - Cosuj 2 




k, 


I 


d 

Ife- 


■ 


Using the results (6.3.4) and (6.3.6) in equation (6.4.27) we get 



E(y)r 


h'(u)du 


H M. 

fLjA 2 ) 


\2a 


Q(e)sin0de 


(Cosy - Cosu) 2 dy b (CosQ-Cosy)^ 


1A 

k dy 


Sin0d0 


(Cos0 - Cosy )2 


i; 


' b E(y)dy 


t£- 


h'(u)du 


(Cosy - Cos0) 2 (Cosy - Cosu) 2 


+/; 


E(y)dy 


g'(u)du 


T + J( 


■e E(y)dy 


'o 1 1 Jo 1 

(Cosy - Cos0) 2 (Cosy - Cosu) 2 (Cosy - Cos0) 2 


c lc 


k'(u)du 


(Cosy - Cosu) 2 


b <0 < c (6.4.28) 


Changing the order of integration in above equation, we get 


E(y>r- 


h'(u)du 


' (Cosy - Cosu) 2 


Ql(y)-i 

71 


fE(y)dy 

L (Cosy-Cosu)^ 


+ 


£ 


■a g'(u)du f a 


(Cost - Cosu) 2 


{ E (y)dy + | 


* k'(u)du 


(Cost - Cosu) 2 


r fE(y)dy 


d ry 


Sin0d0 


dy [(Cos0 - CosyXCost - Cos0)]2 


b < 0 < c (6.4.29) 


1 


where 


J Sin — 


n 2a 


Q(0)Sin0d0 


dy (Cos0 - Cosy) 2 


(6.4.30) 
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d fy SmQdQ __ Siny(Cost“Cosb )2 

dy Jb , l i i (6.4.31) 

(CosG - Cosy) 2 (Cost - Cos0) 2 (Cosb - Cosy ) 2 (Cost - Cosy) 2 


Siny(Cost-Cosb) 2 


in equation (6.4.29), we get 


E(y)j: — -- ' (u)du T =Qi(y)-- — 


(Cosy - Cosu) 2 


(Cosb - Cosy) 2 


fb E(t)(Cost - Cosb) 2 dt rc h'(u)du 

° (C0St “ C0Sy) b (Cost - Cosu)i 

1 

- f aE (tX Cost ~Cosb) 2 dt ra g'(u)du 
° ( Cost - Cosy ) ‘ (Cost - Cosu)i 

r QE(tXCost - Cosb) 2 dt k'(u)du 

“ (C0St - C ° Sy) “(Cost-Cosu)!. 


b < y < c (6.4.32) 


Putting the value of last integral of second term on the right hand side 
of eq. (6.4.32) from eq. (6.4.18), we get 


E(y)f 


h'(u)du 

(Cosb - Cosy) 2 


Qi(y) — 


(Cosb-Cosy) 2 


rb E(tXCost - Cosb) 2 ^ rc h'(u)du 
■*° (Cost -Cosy) ^ b ' 


(Cost - Cosu) 2 


+ f aE(tXCost-Cosb) 2 dt J ^(t) _ rc h'(u)du 

0 -Cosy) W^Cost-Cosu^ 


r^ k f (u)du + r 6E(tXCost - 
^(Cost-Cosu ^ 0 (C0St - 


X 

tXCost - Cosb) 2 
(Cost - Cosy) 




jft k f (u)du 
(Cost - Cosu);> 


b < y < c (6.4.33) 


i 

(Cosy - Cosu )2 


Qi(y)+P 3 (y)- 


_1 Siny 

K (Cosb - Cosy )2 


rb E(tXCost - Cosb )2 dt j-c h'( u)du 
a (Cost - Cosy) •’b 


(Cost - Cosu )2 


1 

j^E(tXCost - Cosb )2 dt rz k'(u)du 
r <“-<*) l< ’ < ‘ 


where 


p ( y ) = sin Y r a fi(t)(Cost - Cosb) 2 dt 
(Cosb-Cosy)j (Cost -Cosy) 


(6.4.35) 


Let B(u, y) = Sin y fb E^XCost-Cosb^dt 
(Cosb-CosyL* (Cost -Cosy) 


(6.4.36) 


and C(u,y) = 


Ja (r, 


1 

ost-Cosb) 2 dt 


(Cosb-Cosy)5 ( Cost ~ Cos y) 


(6.4.36) 


Now equation (6.4.34) can be reduced to the following form 


E(y)^ — '( u)iu i = Qi(y)+Ps(y) 

(Cosy - Cosu )2 





b < y < c (6.4.38) 


Again applying the results (6.3.5) and (6.3.7), in equation (6.4.38) we 


1 d re {Qj (y) + P 3 (y)}Sinsds 1 d 


Sinsds 


E(y)h'(u) 


b < y < e (6.4.39) 


Sinsds 


‘osu - Coss )2 (Coss - Cosv )2 


(6.4.40) 


osc - Coss )2 (Cosv - Coss) 


equation (6.4.39) can be rewritten as 


E(y)h'(u)= Q 2 (u) 


b < u < c (6.4.41) 


where. 


1 d rc {Q 1 (y)+P 3 (y)}Sinsds 
7tdu*’ u a, _ ~ _vi 


(6.4.42) 
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B(s,v) ; 


Sinu(Cosv - Cosc)2 

— _ I 

(Cose - Coss)2 (Cosv - Coss) 


B(v,y) 


(6.4.43) 


and cM=-?=feM_c( Vi ,) 
(Cose - Coss)2 (Cosv - Coss) 


(6.4.44) 


Again using summation result (6.3.3) in terms of integral in equation 


(6.4.10) we get 


fk'(u) + fk'(u) p (u ’ e) m d Udy 

Jd Je ; J Jo [(Cosy-CosuXCosy-Cos0)f' /2 y 


dSin- R(0)- £g'(u)du +£h'(u)du 


rmin(u,0) E(y)dy 

0 (Cosy - CosuXCosy - Cos©) 1 


d < 0 < c (6.4.45) 


J^k'(u)duJ“- _ 

(Cosy - Cosu)2 (Cosy - Cos0)2 

+ £*(»)< ^ r 

(Cosy - Cosu)2 (Cosy - Cos0)2 


re pin- .R(0)- Jg'(u)du+J h'(u)du 


j 6 E(y)dy 

(Cosy - Cosu)2 (Cosy - Cos0);> 


d < 0 < n (6.4.46) 
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Changing the order of integration in eq. (6.4.46), we obtain 


f 9 E(y)dy r% k'(u)du 

id 1 Jy 1 ~ 71 

(Cosy - Cos0 )2 (Cosy - Cosu )2 


r 0\ 2a 


Sin- 


2 ) 


,R(9) 


E(y)dy g'(u)du f t E(y)dy 


r a £tyjgy r a gyujuu 

Jo . . i Jo , . i + Jo 


l0 / v 1 JO , 1 ' JQ . ,1 

(Cosy - CosG )i (Cosy - Cosu)2 (Cosy - CosOji 


cf 

Jb 


h'(u)du 


T + r 


E(y)dy 


rj;- 


h'(u)du 


(Cosy - Cosu )2 (Cosy-Cos0)2 ^ (Cosy-Cosu )2 


+ 


f 

Jo 


E(y)dy 


k'(u)du 


(Cosy - Cos 0)2 (Cosy - Cosu )2 J 


d < 0 < 7i (6.4.47) 


Using the results (6.3.4) and 5.3.6), we get 


N 2a 


E(y){* 


k'(u)du 


d y [sin|J ,R(0)Sin0d0 

“Id r ~ 


(Cosy - Cosu )2 ^ (Cos0 - Cosy )2 


■a E(y)dy ra g'(u)du 


_r 

1 Jv 


+ 


1 d jv Sin0d0 J 
71 ^ (Cos0 - Cosy)-? (Cosy - Cos©);? V (Cosy - Cosu )2 
•b E(y)dy fc h'(u)du E(y)dy 


f 

Jo 


•r- 




(Cosy - Cos0)2 (Cosy - Cosu )2 (Cosy - Cos0)2 


€ 


h'(u)du 


• + 


E(y)dy 


rJT- 


k'(u)du 


(Cosy - Cosu )2 (Cosy - Cos0)2 (Cosy - Cosu )2 


d < y < 7t (6.4.48) 
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Changing the order of integration and applying the result. 


Sin0d0 


Siny(Cost - Cosd )2 


1 I } ]_ 

(Cos0 - Cosy)! (Cost - Cos0)! (Cosd - Cosy)! (Cost - Cosy)! 


(6.4.49) 


in equation (6.4.48), we get 


l.»f k '' UMl: , =R,(y)-l -- Sin >- , 
' (Cosy - Cosu) 2 % (Cosd - Cosy)2 

x f a E ( t )( CoSt ~ Cosd)! dt j-a g'(u)du 
*^° (Cost - Cosy) « 1 


(Cosy - Cosu )2 


rb E(t)(Cost - Cosd)! dt re h'(u)du 

. - T J b “ T 

(Cost - Cosy )2 (Cost - Cosu )2 

l 

jc E(t)(Cost - Cosd)!dt h'(u)du 
(Cost - Cosy)! (Cost - Cosu)! 


fd E( t )(Cost - Cosd) 2 dt pc k'(u)du 

Jo 11 Jd 1 

(Cosy- Cosy ) 2 (Cost - Cosu) 2 


d < y < tc (6.4.50) 


where 


( oy« 

Sin - .R(0)Sin0d0 

R l(y) = ±PL_?i 

d * Jd (Cose-Cosy i 


(6.4.51) 


Now putting the values of last integrals of first and third terms on the 
right hand side from equations (6.4.18) and (6.4.38), we get 
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PjfoXCost - Cosd)2 
(Cost - Cosy) 


(Qt (t)+ P 3 (t)XCost - Cosd )2 dt 
(Cost - Cosy) 


E(tXCost - Cosd )2 
(Cost - Cosy) 


Again applying the results (6.3.5) and (6.3.7) in equation (6.4.52), we 


rbE(tXCost-Cosd)5dt t | 

rcE(t)c(v,t)XCost - Cosd )2 dt 

1 

A .'Ilf 

(Cost - Cosy) J 

b (Cost - Cosy) 

. » . mS , 

j 


E(y)h'(u) = r{ R i(y) +R 2 (y)}Sinsds 
7tdu^ u - - 1 


1 d rx Sinsds 


(Cosu - Coss )2 


+ —j — f 

7C 2 du^ u 


fc U(v,y) h'(v)dv _ fic V(v,y)k'(v)dv 

/ . ^ , J] b l + id I . 

(Cosu - COSSJ 2 [_ (Coss - Cosv )2 (Coss - Cosv )2 J 


(6.4.56) 


Changing the order of integration and using the result 


Sinsds 


Sinu(Cosv - Cos7t)2 


d ju 

(Cosu - Coss )2 (Coss - CoSv )2 (CoS7I - Coss )2 (Cosv - Coss) 


in equation (6.4.56), we get 


E(y)k'(u) = R 3 (u)+ i£h'(v)U(s,v)dv + ±£k'(v)V(s, v)dv 


(6.4.57) 


jr 


d < u < re (6.4.58) 


where. 


r ( u ) = _I A f TC ( R i (y) + R 2 (y )}sm S ds 

7idu^ u /-_ _ \I 

(Cosu - Coss )2 


(6.4.59) 


U(s,v)= Si "“(C°sv-Cos g )2 u(v>y) 
(COSTC - Coss )2 (Cosv - Coss) 


(6.4.60) 


■ ¥» \ 


and 


V(s,v): 


1 

Sinu(Cosv - Cos7t)2 
(Costt - Coss);* (Cosv - Coss) 


-V(v,y) 


(6.4.61) 
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r ; 


Equations (6.4.22), (6.4.41) and (6.4.58) are Fredholm integral 
equations of the second kind which determine g'(u), ff(u) k'(u) respectively. 
Knowing the values of g'(u), h'(u) and k'(u) the values of coefficients A n , can 


be obtained from (6.4.4). 


PARTICULAR CASES 


If we let d-7t in equations (6.2.1) to (6.2.5), they reduce to the 
corresponding quadruple series equations considered earlier in [66] and the 
above solution agrees with that obtained previously. 

Similarly, the solutions of corresponding triple and dual series can be 
obtained as particular cases. 
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TWO GRIFFITH-CRACKS AT THE INTERFACE OPENED 


BY FORCES AT CRACK SURFACES 


7.1 INTRODUCTION 


In the present chapter we are going to find the effect over physical 


quantities (important in fracture design criterion) due to the presence of two 


cracks in dissimilar media 


The cracks occupy the regions y = 0, b<|x|<c at the common interface 


and mathematically the problem is reduced to the following boundary value 


problem, see figures 7.1 and 7.2 


with, b<|x|<c and (±) sign over 0 refer to the quantities corresponding to y>0 


and y<0. There are the following continuity conditions 




Jj^l ll 


Isotropic Medium 


Fig. 7.2: Geometry with Boundary and Continuty Condition 




x-axis 


(-b,0) 


Orthotropic Medium 


tfyy(x,0‘ ) = Oyy(x,0') 


(7.1.6) 


for b<|x|<c. 

The following assumptions are due to the symmetry in crack geometry: 


CTxy(x,(T) = O sy (x,0*), 

b<x<c 

(7.1.7) 

a y y(x,0') - Gyy(x,0'), 

b<x<c 

(7.1.8) 

Uy(x,0 ) = U y (x,0*), 

b<x<c 

(7.1.9) 

u s (x,O f ) = u x (x,0'). 

b<x<c 

(7.1.10) 


We concentrate analysis over 1st and IVth quadrant only. We checked through 

out that 

u y (x,(T)>0 (7.1.11) 

Uy(x,0')>0 (7.1.12) 

which means that the cracks really open out. 

The plan of the chapter is as follows: In section 7.2 the boundary value 
problem will be reduced to triple integral equations. In section 7.3, the 
solution of triple equations will be given in section 7.4, the approximate 
solution of coupled Fredholm integral equations of second kind will be given, 
in section 7.5, the physical quantities in the vicinity of crack tip will be given. 
The surface forces will be assumed as equal, constant and uniform. 
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7.2 REDUCTION TO TRIPLE INTEGRAL EQUATIONS 

The component of stress and of displacement are assumed by means of 
the following equations (in the absence of body force). 

For isotropic case (y>0) 


03 


xx , 


3 X d y 


(7.2.1) 


do v 


*y , ^yy __ q 


d x dy 


(7.2.2) 


yK ,y) ~ tie rw ~ 


" H fe y)-«-(n- 2^ 2 y) 

5y <9y 


d£ (7.2.3) 


u x (x, y )=2(LL!L)r££fe) 


TtE i § 


,a 2 H 


(l-n)— - + T^ 2 Hfe,y) 

oy~ 




(7.2.4) 


with 


Hfey)=[Afe)+yBfe)]e-» 


(7.2.5) 


where A and B are two arbitrary constants. 


For orthotropic case (y<0) 


u y( x >y)= |^|^[aiiG,yyy-^ 2 (a 12 + a 66 )cG,y]^ 
0 s 


(7.2.6) 
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(7.2.7) 


u x ( x >y) = J^~^[ a nG,yy - a 12 ^ 2 G^ 

0 ^ 

Gfe. y) = — [fc - r 2 )Cfe)- D(5)Je r * ] 

r I _ r 2 


(7.2.8) 


where C and D are two arbitrary constants and x\ an d r 2 are two roots of 


r 4 + B,r 2 + B 2 =0 


(7.2.9) 


where 


B, 


2( a | 2 + a 66) g _ 3-22 


a ll 


a ll 


(7.2.10) 


* 

s 

ll 

'I'll 






and an - a 66 etc. are elastic constants of orthotropic medium. The stress- 

strain relations are: 


For isotropic case (y>0) 

2^ =(l-nKx-ntfyyl 

2{iOyy — (l — Tj)o«^ — tIOxx 
2M- e xv =<r xv 


(7.2.11) 


where p. and r\ are modulus of rigidity and poisson’s ratio of isotropic 
medium. 


1 
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For orthotropic case (y<0) 


dux 

- anCJxx +a 12 ayy 


duy 


dy 


a 12 a xx + a 22 a 


yy 


dux duy _ 

+ — ~ a 66 a xy 


dy dx 


(7.2.12) 


The boundary conditions (7.2. l)-7.2.2) and the continuity conditions 
(7.2.7)-(7.2.8), and then getting Fourier inversion of the functions given 


5hcfe)-Dfe)]+Bfe)-5A(5)=0 


(7.2.13) 


Afe)^Cfe)=P,fe)/5 : 


(7.2.14) 


with 


p i fe) = - J (ft (x) - p 2 (x))cos(^x)dx 

71 


(7.2.15) 


Now using the relations (7.2.3)-(7.2.8) and (7. 1 ,9)-(7. 1.10) we get the 
following integral relations: 


J^sin(^xXK 5 C(^)+K 6 D(^)]d^=P 4 (x) c<x<oo, 0<x<b (7.2.16) 


J^cos(^xXK 7 C(^)+K 8 D(^)]d^ = P 5 (x) c < x < oo, 0 < x < b (7.2.17) 
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1 




| 

i 

; j 


(7.2.18) 


K 5 =K 0 -K 1 +2K 0 (l-nXl-r 1 ) 
K 6 =K 2 -2K 0 (l-n) 
K 7 =K () {l-2(l-2T 1 Xl-r 1 )}+K 3 
K 8 =2K 0 (1-2 ti)+K 4 

v — MLiHl j< — r~a —a 
7CE 

K 2 = r 2 a i i > K 3 = rj(a u rf -a 12 -a 66 
K 4 = a ll r 2 _a 12 ~ a 66 


„ | Pife)sinfex)d^ 

l $ 

p Jf Pife)cosfex)d^ 

5 J ? 

o ^ 

Similarly we now evaluate 

axy(x,0 + ) = axy(x,0") = 0, b<x<c 

ayy(x,0 + ) = ayy(x,0‘) = -(Pi(x) + P 2 (x)), b<x<c 

Thus above two relations along with (7.2.3)-(7.2.8) and the 
corresponding stress-strain relations lead us to 

00 

fe 2 sm(^t)[cfe)]dfe)= P 6 to b < X < c (7.2.21) 


(7.2.19) 


(7.2.20) 



jV cosfexjcfejdfe) = P 7 (x) b < x < c (7.2.22) 

0 


with 


1 00 

p e( x ) = I“^r Jp^sin^x^x (7.2.23) 

p 7 (x) = -- Jpi £ cos(§x)+ - (pj (x)+p 2 (x)) (7.2.24) 

2 o 2 

Thus the boundary value problem is reduced to the solution of triple 
integral equations (7.2.16)-(7.2.24). 

7.3 SOLUTION OF TRIPLE INTEGRAL EQUAITONS 
We assure that 

5[K 5 Cfe)+K 6 Dfe)] = ^) (7.3.1) 

5[K 7 Cfe)+K 8 Dfe)]= ¥ fe) (7.3.2) 

where (j> and \) / are two new functions 

Then the integral equations (7.2.16)-(7.2.17); (7.2.21)-(7.2.22) become 
as new triple integral equations. 

OD 

J<l>fe)sinfef)d£ = P 4 (xi 


x e Ij u I 3 


(7.3.3) 



J v(£)cos(£x)d£ = P 5 (4 x s Ij u I 3 (7.3.4) 

0 


and 


J^[ K i 2 <l>(^)+K : 1 3 \[/(^)]sm(^x)d^ = P 6 (x), x e I 2 (7.3.5) 

0 


j £[ K 1 o<KO - K i 1 v|/(^)]cos(^x)d^ = P 7 (x), x € I 2 (7.3.6) 

0 


where K !0 and K 13 are given by (3.3.7) and 


Ii = [0,b], h = [b,c], I 3 = [c,oo] 


(7.3.7) 


Now to solve the system of equations (7.3.3)-(7.3.6) we make use of 
the method of Srivastava and Lowengrub [52], We assume that 



(7.3.8) 


%( 


54© =2 


(b 


| h(t)sin £tdt- J+ J P 5 l (t)sin £t dt 


VO c J 


(7.3.9) 


Then the equations (7.3.3)-(7.3.4) are satisfied identically if 


| g(t) = P 4 (b) - P 4 ( c ) 


(7.3.10) 




(7.3.11) 


Substitution of <j>(£) and y(£) from (7.3.8)-(7.3.9) into (7.3.5)-(7.3.6) 
and evaluating the integrals we get 


(7.3.12) 


(7.3.13) 


where 


(7.3.14) 


(7.3.15) 


Now rewritting (7.3. 12)-(7.3.13) into the following form as 


(7.3.16) 


PioW=P 8 W + K 10 g(x) 
Pn(x)=P 9 (x)-K 13 h(x) 


(7.3.17) 


Inverting (7.3.16) by using Srivastava and Lowengrub method [52] 


h(t): 


7tK n 8(t) 


V 

f 


xS(x)P 10 (x)dx 

2 .2 
X -t 


+ Dl 


(7.3.18) 


8(0 = 


TClC ] 28 (t) 


fx8(x)P n (x)dx 

2 "2 + L >2 

b X -t 


J 


(7.3.19) 


^(t)= ft 2 -b 2 )(c 2 -t 2 )] /2 


(7.3.20) 


where Di and D 2 and two arbitrary constants which will obtained through 
(7.3.1 0)-(7 .3.1 1) and (7.3. 1 8)-(7.3. 1 9). Substituting the values of Pio(x) and 
Pn(x) from (7.3.17) into (7.3. 1 8)-(7. 3.19) and evaluating the integrals. 


h(t)= P] 2 (t)+ 


Kiq | x8(x )g(x)dx 
x 2 -t 2 ' 


(7.3.21) 


g(t)= Pi s(t) " 


K 13 t < r 8(x)h(x)dx 
7tKi 2 8(t)| x 2 -t 2 


(7.3.22) 


Pl 2 (t) : 


D,+| 


xS(x)P 8 (x)dx 

2 72 
x - 1 


7tK n 8(t) 


(7.3.23) 


Pi 3 (t)= 


D, + f 


r S(x^P 9 (x)dx 


7 ' i „2 +2 


1 


•rrTC^Sft) 


(7.3.24) 





7.4 REDUCTION TO AND SOLUTION OF FREDHOLM INTEGRAL 
EQUATIONS 

We substitute for g(t) from (7.3.22), (7.3.24) into (7.3.20) and evaluate 
certain integrals, we get 


h (t) = Pi 4 (t)- ^ J^»?H(t,a)da 
5 W i or - t 2 


+ llLr 
2S(t ){ 


S(a)P 9 (a) 

2 .2 
a -t 


H (t,a)da 


with 


H 0 (t)=itlog 


c-t b + t 


c + t b-t 


(7.4.1) 


(7.4.2) 


(7.4.3) 







and 


H(t,a) = H 0 (t)-H 0 (a) 


K m 


27i 2 K H K 12 


K 


K 


15 


10 


2% K„K 12 


Similarly 



(7.4.4) 
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*,(t)H5(a) 


(7.4.6) 


(7.4.7) 


(7.4.8) 


(7.4.9) 


Thus knowing the solution of (7.4.1), we will substitute in (7.4.6) to 
get the solution of g(t). 


We assume h(t) as 


(7.4.10) 


Substituting this value of h(t) in (7.4.1) and then comparing the coefficient of 
{6 n } we get 


H O (0=P M (t) 

H l(t)=-K, 5 I„(t) 

^2(0= - Ki 5 li(t)+r|Ki 5 I 0 (t) 
H 3 (t)=-K, 5 l2(t)+r 1 Ki5(l + Ki5)li( t )-r ] Ki 5 I 0 (t) 


(7.4.11) 


(7.4.12) 


(0 JKq (a, t) I n _| (a)da 


b 


I 0 (t) = jK 0 (a,t)P 14 da 

b 


Thus h(t) is written as 



h(t) = P, 4 (t) - 5K ! 5<|) 2 (t) 


with 


^2 (t) - K 20 Io (t) + K 2 1 1] H + K 22 12 (t ) 


and 


K 2 o - l-rjd-rjd" 

K 21 =8K 15 -5 2 ri (l + K 15 ) 
K 22 =8 2 K 15 


where I 0 , 1 1? I 2 etc. are given by (7.4.12)-(7.4.13). 

Substituting this value of h(t) from (7.4.15)-(7.4. 16) into 
evaluating certain integrals, we get 


(7.4.13) 

(7.4.14) 

(7.4.15) 


(7.4.16) 

(7.4.6) and 



where 


. . ... .., 




<J> 3( t )=4>i(t)[Gf ) (t)-§K 15 K 20 Gf )) (t)+8K 15 K2 1 G} 0 )(t) 
+ SK 15 K 22 0|’)(t)]-G< , l(t)-8K l5 K 21 G 1 W(t) 

+ 8K, 5 K 20 G^(t)-8K I5 K 22 G^(t) 


(7.4.17) 


and 


GL"»(.)=J^Tl m -M 2 ”da 

a J 


,a~ -t 


I 3 (t)=P,4(t) 


(7.4.18) 


while 


<j>i(t) is given by (7.4.9). 


and the arbitrary constants D] and D 2 are evaluated through (7.4. 15)-(7.4. 18) 

and (7.3.1 0)-(7.3. 1 1) 


D] =rcK n [P s (b)-P 5 (c)-J 3 3+J 4 ]/J 


31 


^i!^[P 4 (b)-P 4 (c)-K, 6 J 5 ] 

J 32 


(7.4.19) 


D, 2 = 24 P 4 (b)-P 4 (o)-Ki 6 J 5 ]K 


12 


(7.4.20) 


1 

J 3 i = -F(7c/2,Pi ) = Complete elliptic integral of first kind 
b 


K( 


(c 2 -b 2 )b : 


J 






7 i 
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n is elliptic integral of third type 


7.5 PHYSICAL QUANTITIES 


The quantities, which are important in fracture design criterions, are 
the components of stress in the neighbourhood of crack tips and the crack 
opening displacement. The components of stress are evaluated through the 
following relations. 
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Stress Components 


I 

r 

; Shear Stress 

a xy (x,0 + )-a xy (x,(r)=0, xeljulj (7.5.1) 

°xy ( x ’0 + )+ ^xv (x,0“ ) = (1 +r 1 )Jsin^x^ 2 Cfe)-P 1 fe)]d^ (7.5.2) 

o 

The above integral becomes as 

<?xv ( x > 0+ )+ Oxv ( x >°~ )=( 1 + I i)J [ K i 2<i>fe) + 2 vfe)l 

o (7.5.3) 

sin(^x)d^ - P 6 (x), x e Ij u I 3 

Using relations (7.3.8)-(7.3.9) and then changing the order of 
integration and evaluating the integrals. 


a xy (x,0 + )+a xy (x,0")= K 12 (l + ~ J + J 

L b l ~ x n vo c ; 


P 4 (t)dt 


-(l + r 1 )K 13 P 5 (x)-P 6 (x), xel x KjI 3 (7.5.3) 

Normal Stress 

Similarly the normal component of stresses are given as 
a xy (x,O + )- 0 xy (x,O“)=O, xeli^l 3 (7-5.5) 

and 
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Now substituting for g(t) and h(t) from (7.4.15)-(7.4.18) into (7.5.4) 
(7.5.6) and evaluating the integrals, we get 


2 x (P 15 (x)+ K 16 0 3 (x)) 


{P 14 ( x )-6K, 5 4>,(x)}+f 0 (x) 


(7.5.6) 


F„(x) = K,„P 4 (x)+- J+J P‘(t) 


,2 2 
t -X 


%lwl» 


t -X' 


•P 7 (x) xel, ul 3 


Displacement Component 


The displacement components are evaluated through relations 
(3.2.3 );(7. 2.2) - (7.2.3);( 7.3 . 1 ) - (7.3.2);(7.3.8) - (7.3.9) and evaluating the 
integrals we get 


iso u y (x,0)= u y (x,0 + ): 


>(l + r\) 2(l + t))P 5 (x) 


• + r/c 


g(t)- J+J |pi(t) >log|x 2 - 1 2 dt 


+ K 1 + , Jh(t)dt-P 4 (c) 


,s + tr r l K 8 _K 7“ r l K 6 
^10 _ k (j - 


Kjj — K 0 


(l-2n)K 5 -r 1 K 6 


(7.5.8) 


Further 


orth. u y (x,o) = u y (x,0 + )= Jcosfex)^(K 3 C(^)+ K 4 Dfe))]d£ 


= J C0S (^ X ( K 3 ^ - ~ ( K 7* “ ] K 5 V) 

0 L K 9 r^9 


K c 


Thus 


u y (x 3 °-)= Jcosfex^ +K 1 + 4 M/]d^ 


K 


+ _k 3 k 8 -k 4 k 7 


10 


Kc 


K IT' \S 1 / 

„ 4 K 5 -k 3 k 6 

14 K„ 



u y ( x >°~)=Ps(x)+— 

71 


Jg(t)logx 2 -t 2 


dt 


00 ^ 

J+j pi(t)io g 

U c 


2 .2 
X -t 


dt 


(7.5.11) 


+ K, + ! 


Jh(t)dt + P 5 (c) 


xeh 
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Similarly 



; (x,0 + )=P 4 ( x )K 0 +Kji fg(t)dt + P 4 (c) 


• + C f c co 

b vb c J 


Kj^Kg + 2(1 - r|)K 2 


K, + 5 K 6 + 2(1 - r|)K 


*17 = K 0 


*15 = 1 + 2 (l - rjXl - r, ) 


u x (x,0“)= fecosftxl Cfe)- a ^ D d£ + K 0 P 4 (x) 


K 8 <j)-K 6 \|/\ K 7 ^-K s vj/ 


jsin(£xj K 


J sin^x^K^)- K, + 9 vj/(^)]d^ + K 0 P 4 (x), 



M x ’° j=KoP 4 (x)+Kf 6 Jg(tjtit + P 4 (c) 


+ - 


Kj + 9 


K 


c f b oo A 

|h(t)- J+j Pj(t)>log| 
b VO c J 


(7.5.15) 


x-t 


x + t 


dt 


X € I« 


7.6 SPECIAL CASE 


The given forces at crack faces are assumed to be constant, uniform 
and equal, see figure 7.3 

Pi(x) = P 2 (x) = P 0 (7.6.1) 

then 


Pi(5) = 0 


and 


P 4 (x) = P 5 (x) = P 6 (x) = 0 


P 7 (x) = P 0 


( 7 . 6 . 2 ) 


(7.6.3) 

(7.6.4) 


then (7.3.10)-(7.3.1 1) reduce 




Pio(x) = Kiog(x) + P 0 


(7.6.8) 


Pn(x) = -Knh(t) 


P 12 (t)- 


Dj +p 0 ^t 2 -c 2 -b 2 ^ 


JtK n 8(t) 


P, 3 (t)- tD2 




P 14 (0 - p 12 (0+ 2 {c - b + H 0 (t)} 


P| 5 (t) 


7iK 12 8(t) 

HoW 


D- 


K 


13 


7CK 


11 


( Dl + ^(c 2 -b 2 )} 


2t" 


+ 7ip 0 ( c -b + H 0 (t) 


where (7.6.3)-(7.6.7) are substituted in (7.4.19)-4.4.21), we get 


Dj = ^li[j 4 - J 33 ]+2 tc 2 K 11 K 12 - -15 . K 16 J 32 J 5 


'll 


'31 


J 33 =-l^L-2b 2 J 31 ] 
33 7cK n 1 31J 


(7.6.9) 


(7.6.10) 


(7.6.11) 


(7.6.12) 


(7.6.13) 


(7.6.14) 


The constants J 4 , J 5 , J 32 , J 31 can be evaluated numerically, through 
these integrals which are to be handled carefully for numerical integrations. 


LS'l 


{/.o. id ; 



The constant J32 has logarithmic oscillations, hence the evaluation has 


to be done very carefully. The functions h(t) and g(t) given by (7.4.15) 


(7.4. 16) with P 1 4 , P 1 5 given by (7.6.12), (7.6.13) and constants Dj and D 2 the 


there in by (7.6.1 4)-(7. 6. 15) and remaining constants and functions as defined 


previously 


Thus the physical quantities can be evaluated by numerical integration 
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